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A superposition of a matrix ensemble refers to the ensemble constructed from two indepen¬ 
dent copies of the original, while a decimation refers to the formation of a new ensemble by 
observing only every second eigenvalue. In the cases of the classical matrix ensembles with 
orthogonal symmetry, it is known that forming superpositions and decimations gives rise 
to classical matrix ensembles with unitary and symplectic symmetry. The basic identities 
expressing these facts can be extended to include a parameter, which in turn provides us 
with probability density functions which we take as the definition of special parameter de¬ 
pendent matrix ensembles. The parameter dependent ensembles relating to superpositions 
interpolate between superimposed orthogonal ensembles and a unitary ensemble, while the 
parameter dependent ensembles relating to decimations interpolate between an orthogonal 
ensemble with an even number of eigenvalues and a symplectic ensemble of half the number 
of eigenvalues. By the construction of new families of biorthogonal and skew orthogonal 
polynomials, we are able to compute the corresponding correlation functions, both in the 
finite system and in various scaled limits. Specializing back to the cases of orthogonal and 
symplectic symmetry, we find that our results imply different functional forms to those known 
previously. 


1 Introduction 

Dyson ^ introduced three ensembles of random unitary matrices — the circular orthogonal ensemble 
(COE), circular unitary ensemble (CUE) and circular symplectic ensemble (CSE). The corresponding 
joint eigenvalue probability density functions (PDFs) were calculated to be 

2 n (1-1) 

l<j<k<n 

where C is the normalization (throughout the symbol C will be used to denote some normalization) and 
/3 = 1 for the COE, /3 = 2 for the CUE and /3 = 4 for the CSE. In the theory of the COE, a technique 
known as integration over alternate variables has a special place. This technique draws one naturally 
to study statistical properties of every second eigenvalue (parity respecting correlations) as well as the 
statistical properties of the complete COE sequence (parity blind correlations). Thus for matrices from 
the ensemble COE 2 n (2n x 2n members of the COE) Mehta and Dyson considered the statistical 
properties of every second eigenvalue (referred to as an alternating sequence) by integrating out the 
complementary alternating sequence. With the resulting distribution denoted alt(COE 2 „) they showed 

alt(COE2„) = CSE„, (1.2) 

where on the RHS we mean the joint eigenvalue distribution of the ensemble CSE„ (below we will state 
similar equations with this convention without further comment). 
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Prior to the work of Mehta and Dyson, Dyson was led to conjecture that superimposing two 
independent eigenvalue sequences from the COE and integrating out every second eigenvalue leaves an 
eigenvalue sequence with the same joint distribution as the CUE, and thus 


alt(COE„ U COE„) = CUEn- 


(1.3) 


This was subsequently proved by Gunson The results o and (OD together imply that the 

physically important gap probability, that is the probability that the interval [—s, s] is free of eigenvalues, 
is inter-related for the three ensembles COE, CUE and CSE. As formulas for this quantity were known 
for the COE and CUE in terms of the Fredholm determinant of the integral operator on [—s, s] with 
kernel sin7r(a: — y)/'K(x — y), the inter-relationships imply a formula for the gap probability in the CSE 
in terms of the same Fredholm determinant . 

Some years after the pioneering work of Dyson and Mehta, Baik and Rains 0 were led to study 
the distribution of every second row of random Young tableaux specified according to some specific 
probability measures. In particular, a probability measure was identified for which the distribution of the 
even numbered rows is independent of a parameter a occurring in the measure. In the continuum limit, 
in which the integer valued row lengths go over to continuous valued variables, the PDF of the 2n row 
lengths is specified by 

^ 2n n 

n (i-4) 

j—l j=l l<j<k<2n 


where 

Xi > X2 > ■ ■ ■ > X2n > 0, (1.5) 

C is the normalization and A is the analogue of the parameter a. For (0) to be normalizable we must 
have A < 1. In the case A = 0 this coincides with the joint distribution of the eigenvalues in the matrix 
ensemble LOE — the Laguerre orthogonal ensemble with parameter a = 0 (for general parameter a the 
Laguerre weight is a;“e“^/^). We recall the LOE distribution of 2n variables with parameter a = 0 is 
realized by the eigenvalues of Wishart matrices X^X where the real matrix X has dimension (2n-|-l) x 2n 
and independent elements, identically distributed with the standard normal distribution N[0,1]. Let us 


denote the ensemble corresponding to (1^) by LOE^„. With ME denoting a general matrix ensemble 
and the operation even(ME) denoting the distribution of the even labelled coordinates with the ordering 
(0, and thus the operation of integrating out the odd labelled coordinates, the result of Baik and Rains 
gives 

even(LOE^„) = LSE„ (1.6) 


where the PDF for the LSE„ is the A ^ —oo limit of the PDF (1.4), and is thus given by 


n 


1=1 


l<j<k<n 


(1.7) 


after re-labelling coordinates {x 2 j-i,X 2 j} > Xj (j = l,...,n) (see ( 3.16D ). As a matrix ensemble the 

LSE refers to the Laguerre symplectic ensemble with parameter a = 0. It can be realized as a 2n x 2n 
antisymmetric matrix, in which the elements are pure imaginary numbers with each 2x2 block having 
a real quaternion structure. Such matrices are equivalent to block matrices of the form 


Onxn 

^nxn 


^nxn 

Onxn 
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where X is an antisymmetric complex matrix. As is evident from the transition from (III) to (|h^ , the 
corresponding eigenvalue spectrum is doubly degenerate. 

Implicit in Q is the matrix ensemble corresponding to the eigenvalue PDF 

- 2n n 

ix2j-i-X2k-i){x2k-X2j), (1.8) 


i=i 




l<_ 7 <A;<n 


where the ordering (^) is assumed, and as with (1^) we must have A < 1 for the PDF to be normalizable. 
When A = 0 the matrix ensemble with this eigenvalue PDF is the superimposed ensemble 

LOE„ U LOE„ 


(two independent copies of the LOE). The A ^ —oo limit of (1.8) gives the LUE, parameter a = 0, joint 
eigenvalue distribution 

1 " 

IIP II I ' 

C 


(1.9) 


■""" n (xk-xjf 

j—l ^^j<k<n 

where as in (0) we have re-labelled the coordinates {x 2 j-i,X 2 j} Xj {j = l,...,n). As a matrix 
ensemble the LUE refers to the Laguerre unitary ensemble, which can be realized by matrices of the 
form X^X where X is a n x n matrix with independent, identically distributed complex Gaussian entries. 
Analogous to the derivation of ( [l.4[) , one can use arguments based on the underlying combinatorial model 
to show that 

even((LOE„ U LOE„)^) = LUE„ (1.10) 

independent of the parameter value A. 

Inspired by the results (0 and (|1.10D , Forrester and Rains jl^ considered general matrix ensembles 
with orthogonal symmetry OE„(/) corresponding to the joint distribution 


n (xj-Xk)- 


( 1 . 11 ) 


i=i 


l<j<k<n 


They sought to classify all differentiable weight functions / in ( 1.11 ) such that 

even(OE„(/) U OE„(/)) = UE„( 5 ) (1.12) 

where UE„(g) is the matrix ensemble with unitary symmetry corresponding to the joint distribution 

1 " 

■^Yidixj) n {xj-Xkf. 

j—l l^j<k<n 

Up to linear fractional transformations only two pairs of weights (/, g) were found to possess this property: 


if,9) = 


(e e ^), X > 0 

((1 — (1 — x)°‘), 0 < X < 1. 


(1.13) 


The first of these is the a = 0 Laguerre weight in (lA), while the second is an example of the Jacobi 
weight x^(l — a;)“ with parameter 6 = 0. Furthermore, it was proved |j^ that the statement ( 1.12| ) is 
equivalent to the statement 

even(OE 2 „(/)) = SE„((g/ff) (1.14) 

where SE„(6) denotes the matrix ensemble with symplectic symmetry corresponding to the joint distri¬ 
bution 

1 " 

Y[h{xj) (xj-Xk)^. (1.15) 


C 
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Our interests in this work are parameter dependent generalizations of the orthogonal ensembles, and 
superimposed orthogonal ensembles, specified by the weights / in (1.13). For the weight f{x) = 

(x > 0), these are precisely the ensembles with PDFs ( |l.4D and (1^). For the weight f{x) = — 

0 < a:: < 1, the parameter dependent ensembles are specified by the PDFs 


I zn n 

7= 11(1n( 


i=i 


1=1 


1 - X21-1 \ -^1'^ 
1 - X2l 


l<j<k<2n 


and 


- 2n n ^ 

in(i-xi)<“-‘''*n(T 


1 - X2l-l\-^l‘^ 


1=1 


1=1 


- X2l 


ix2j-l - X2k-l){x2j - X2k) ( 1 - 17 ) 

l<j<k<n 


where, for convenience, we shift the origin and scale the variable x relative to that used in (1.13) so that 


— l<Xj <1 


(j = l,...,2n). 


For ( 1.16| ) and (1.17) to be normalizable we must have A < a + 1. We remark that after shifting the 
origin back to that used in ( 1.13| ) by Xj \{xj + 1), then scaling the variables and parameters 

Xj I—^ XjjL^ Qj I—^ I—^ LA 

and taking the limit L oo, (1.16) and (1.17) reduce to (1.4) and (O) respectively. 


We note that with respect to each even labelled coordinate the parameter A in ( 1.16 ) and ( 1.17 ) 
can be viewed as a change of parameter a a + A in the Jacobi weight. Let us make this change in 
( |1.13| ) and form the identities ( 1.12| ) and (1.14). If we then divide both sides of the resulting identities 
byni:i(i — so that in each case the RHS is independent of A, we deduce that 


even(^OE„(/o, /e) U OE„(/o, /e)) = UE„(g), 

even(^OE2n(/o,/e)) = SE„ (^(g//eU=o 


(1.18) 

(1.19) 


with 


ifo, fe, g) = ((1 - (1 - (1 _ a:)“), -1 < X < 1. (1.20) 

In (p^) , assuming the ordering xi > X 2 > ■ ■ ■ > X 2 n, OE„(/o, fe) U OE„(/o, fe) refers to the ensemble 
with PDF 

n 

Y\_fo{x2j-l)fe{x2j) {X2j - X2k)ix2j-1 - X2k-l), (1-21) 

j=l ^^j<k<n 

while in (1.19|) OE 2 n(/o, fe) refers to the ensemble with PDF 


Y[fo{x 2 j-l)fe{X 2 j) (Xj-Xk). 

j—1 l<j<k<2n 


( 1 . 22 ) 


An analogous argument starting with the Laguerre case of (1.13), modified so that x therein is replaced 
by (1 + A)x, can be used to deduce (0 and ( ll.lOl ). 

The parameter dependent PDFs (|l.4| ), (1.8), ( 1.16 ) and ( 1.17 ) each have at least two interpretations 
in distinct applied settings. One, already made explicit Q, || in the case of (1.4), is as the continuum 
limit of certain measures on partitions. These measures in turn are intimately related to increasing 


subsequence problems 


I , growth models 


and non-intersecting lattice paths 




The other is as the eigenvalue PDF for certain computable ensembles of random matrices. Development 
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of these settings will be undertaken in a separate publication Here we seek the evaluation of the 
multi-point correlation functions, both parity aware and parity blind (recall the terminology from the 
first paragraph) associated with the probability densities (1^), (^), ( 1.16| ) and ( 1.17 ). Let us first revise 
the definition of a multi-point correlation function. Consider a general PDF 


Pixi, 


5 2 /I 5 


5 Un 


(1.23) 


which is symmetric in {xj}j=i^,,,^n and with the support of p on some region I of the real 

line for each coordinate. This can be thought of as the PDF for a two species system of particles free to 
move on I. The (fci, A: 2 )-point correlation function for ki particles of species x and k 2 particles of species 
y is defined as 

P(fei.fc 2 ){a:i> ■,yk 2 ) = n{n -l)---{n-ki + l)n{n - 1) • • • (n - fe -H 1) 

n n „ „ 

^ n n / dxi / dyiip{xi,...,Xn-,yi,-.-,yn)- 

l=ki + ll'=k 2 + l'^^ 

All the PDFs (0), (0), ( |l.l6[) and ( |l.l7| ) can be written in the form (|l.23| ) with xi,.. .,Xn denoting 
the odd labelled particles in the ordering (0)> and yi,... ,yn denoting the even labelled particles in the 
ordering (1.5). 


In the parity aware cases, it follows from general formulas due to Rains that for the PDFs (1.8) 
and (1.16) the (fci, fc 2 )-point correlation functions have the determinant structure 

P(A;i,A;2)(^l7 ■ ■ ■ • ■ ■ 52 /^ 2 ) 

{d^oo{Xj , [d^oe{Xj , y/)] j = 

^ = 1 k2 

[d^eoiyj j Xl)] j=l,...,k2 [Keei^yj , y/)]y ^^3 


= det 


(1.24) 


for certain functions Koo, Koe, Keo and ATee- The latter are defined in terms of arbitrary polynomials 
Pj{y), Qj{x) of degree j, as well as the inverse of the matrix 




dye-‘^^+^'>y/^p,{y) / dxe-^^-^^^/^Qk{x) 

Jy 

in the case of ( |l.4|) , and the inverse of the matrix 

'' dy(l-y)(“-^-')/2^, (y) [\x {1 - x)^‘^+^-^'>/^Qk{x) 






(1.25) 


(1.26) 


-1 Jy 

in the case of (06). Also, a general formula of [Q gives that for the PDFs ( |l.4| ) and ( 1.17 ) the (fci, ^ 2 )- 
point correlation functions have the quaternion determinant (the definition of a quaternion determinant 
is revised in Section 0) structure 


P{k\ ,^2) (^17 ■ ■ ■ j 7 yi 5 ■ ■ ■ 5 yk 2 ) 

[/oo(3^j 5 ^/)]fci [foe {Xj 1 y/)] J = i. 

1 = 1,...,k2 

\feo(yj , X/)] j=i,...,fc 2 [/ee(yj 5 y/)]j ,^2 


= qdet 


(1.27) 


where the /sis 2 are 2x2 matrices with elements defined in terms of arbitrary jth degree polynomials 
Rj{x), together with quantities which differ in their specification depending on whether one is considering 
(0|) or (03). For (y) these quantities are 




(1.28) 
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and the inverse of the antisymmetric matrix 




'-Jo 


- j,fe=0,...,2n—1 


while for (1.17) they are 




(1.29) 


and the inverse of the antisymmetric matrix 

f\l - x)^‘^-^^^^Rjix)<^l{x) dx 




Regarding the parity blind correlations, pk{xi,... ,Xk) (here the symbol Xj is used in its original sense 
of (1^), (1^), ( 1.16 ) and (|l.l7|) ) we note that because (^) and (|l.l7| ) do not vanish when X 2 j-i = X 2 j, 
Pk does not vanish at coincident points and so cannot have a simple determinental form analogous to 
(|1.16). However (|l.4|) and (1.16) can both be written in a form involving a Pfafhan factor, 


2n 


]^w(a;j) {xj - xi)Fi[e{xj,xi]j^i=i^,,,^2n, 

j^l l<j<l<2n 


(1.30) 


for which the general structure of the fc-point correlation is known pi[| in terms of a quaternion determi¬ 
nant 

Pk{xi,. ..,Xk)= qdet[f{xj,xi)]jj^i^,,,^k (1-31) 

where / is a 2 x 2 matrix defined in terms of the same quantities as those specifying (S- 

We find that in the case of the parameter dependent Laguerre ensembles O) and (|l.8D, the corre¬ 
lations can be written in terms of 

n— 1 ^ 

K}:{x,y) :=e-("+^)/2^^L,(x)L,(y), (1.32) 

1=0 

where Li{x) denotes the Laguerre polynomial of degree I and parameter a = 0, with the orthogonality 
property 

POO 

/ e~*Lj{t)Lk{t)dt = hfSj^k, = ^- (1.33) 

Jo 

The function is familiar as determining the correlation function of the Laguerre unitary ensemble with 
parameter a = 0. The latter has the eigenvalue PDF ( |l.9| ). Explicitly, the fc-point correlation function is 
given in terms of (1.32|) by 


Pk{xi,...,Xk)=det K^{xj,xi) 


J j,l=l,...,k 


(1.34) 


Similarly, in the case of the parameter dependent Jacobi ensembles ( 1.16 ) and ( 1.17 ), the correlations 
can be written in terms of 


n—1 -| 

Kilix,y) := (1 - x)(“-i)/^(l - ^ p/“.o)(^)p(“.o)(y) 

1=0 


(1.35) 


where p[°'’^\x) denotes the Jacobi polynomial of degree I and parameter 6=0, with the orthogonality 
property 

/ I oa+1 

(1 _ t)“pj“>'’)(t)p^“>'’)(t)dt = 6/6,-fc, 6/ = , . , _ ^ . (1.36) 


2j + a + 1 
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We remark that 


Pk{xi,. ..,Xk)= det (1 - xi)K:^{xj,xi) 




(1.37) 


is the fc-point correlation function for the eigenvalue PDF 

1 


C 


1<j<fc<n 


Z =1 


which corresponds to the Jacobi unitary ensemble with parameter 6 = 0. 

Let us now turn to the plan of the paper. In Section 2 we take up the problem of the explicit 


computation of the entries of (1.16), and in Section 3 we compute the explicit form of the entries in 
( |l.27 ) and ( 1.31 ). Various scaled limits of the correlations are computed in Section 4, while we finish in 
Section 5 by relating the correlation functions found in this study to correlation and distribution functions 
known from previous studies. 

2 Correlations for superimposed orthogonal ensembles with a 
parameter 

Consider a PDF of the form 


Wo(a;j)we(?/j) {xj-Xk){yj-yk)<iet[K{xj,yk)]j,k=i,...,n, 

j—1 


( 2 . 1 ) 


where Xj,yj G R (although the support of Wo, We may be some subset of R). A general formula of 


gives that the {ki, fc 2 )-point correlation function is given by ( 1.24 ) with 

Tl —1 /.OO 


^OO ^ ) 

Koeix,y) 

Keoiy,x) 

Keeiy,y') 


Wo{x)Qj{x) / k{x', u)we{u)pk{u)du 

j,k=o 

n—l 

Wo{,x)Qj{x) M~* We{y)pk{y) 

j,k=0 

-K(x,y)+ Y { K{u,y)We{u)Qj{u)duj M~f! ( / K{x,v)Wo{v)pk{v)Yv 

\ 

Yi ki{v,y)wo{v)Qjiv)dv] M~,! We{y')pk{y') (2.2) 

: 7 n ' J —OO ' 


j,k=0 


where * denotes the operation of taking the transpose of the inverse, Pj{y) and Qj{x) are as in ( 1.25 ) 
and [Mjk] is the matrix with entries 

/ OO pOO 

dxWo{x)Qj{x) / dyK{x,y)we{y)pk{y)- (2.3) 

-OO j —OO 

This result is relevant to (1.8) and (1.16) because both these PDFs can be written in the form (13) with 


^A{x Laguerrecase 


i^{x,y) = { ,, 

- J Xx>y, Jacobicase 


(2.4) 


where xt = 1 if T is true and 0 otherwise, a fact which can be seen by making note of the following 
determinant identity. 
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Lemma 1. For the orderings 


Xi > ■ ■ ■ > x„ 


yi> ■■■> Vn 


we have 


det[Xxj-y^ >o]i,fc=l,...,ra — Xxi>yi>--->xn>y„- 


(2.5) 


( 2 . 6 ) 


Proof. For the ordering x\ > yi > ■ ■ ■ > Xn > Vn the determinant is triangular with I’s down the 
diagonal, so (|]^) is correct in this case. All other orderings must have at least two a;’s (or two y's) in 
succession. The corresponding rows (or columns) in the determinant will then be equal so the determinant 


vanishes. 

Thus with K{x,y) given by (2.4), we can substitute (2.6) times 

1 — Xi \ -^/2 

i=i t=i 

for the determinant in (|2.lD, and 


□ 


n n - 

'^g^A{xj-yj)/2 Laguerrecase, ^ 

n = ^ yi 


Jacobi case 


Co (a;) = We{x) = 


e (a; > 0), Laguerre case, 

(1 — (—1 < a; < 1), Jacobi case 


(2.7) 


provided the ordering (^) is assumed. On this latter point, (0) is a symmetric function of the a:’s and 
y’s (separately), so the ordering constraint (^) is in fact irrelevant. 

For general PDFs ( |2.1| ) there is of course no explicit formula available for the inverse of [Mjk\ which 
is required in (^). However, for the particular PDFs (|l.8D and ( 1.16 ), and thus K{x,y) given by (2.4), 
this problem can be overcome by choosing {pj{y)} and {Qk{x)} to have the biorthogonality property 


poo 

dywe{y)pj{y) / dxWo{x)Qji{x) = UjSjx- 
> J V 


( 2 . 8 ) 


These polynomials are given simply in terms of the orthogonal polynomials {Lj{x)} in the Laguerre case 
and {P-“’°^(a;)} in the Jacobi case. 


Proposition 2. The sets of polynomials {Lj(a:)}j=o,i,... {Qj"(a:)}j=o,i,... where 

have the biorthogonality property 

poo pOO rt 

J dxe~''^~^'>^^'^Qg{x) = AfpSp^g, ~ 

while the sets of polynomials {Pj'°’°^(a;)}j=o,i,... and {(5^(a;)}j=o,i,... where 

have the biorthogonality property 

J dt (1 - t)<^“+^“^^/^Pp(t) y dx {1 - x)^‘^~'^~^'>^'^Qq{x) = AfpSp^g 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


/-I 

■= 

J \ p . 


2a+l 


(a + 1 + 2p)(p + (a — v4 + l)/2) 


( 2 . 12 ) 
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Proof. Substituting ( ^.9[) in the LHS of (|2.10| ) and making use of (|l.33| ) gives the RHS of ( 2.10 ). Thus 
it only remains to check that Qj{x) is indeed a polynomial of degree j. This latter point follows by 
inspection of (2.£). The verification of (2.12) is done analogously. □ 

The biorthogonality properties ( ^.10 ) and ( 2.12 ) allow in (^) to be replaced by 
The double sums in (2.2) then collapse to single sums. Furthermore, taking note of the values of the 
normalizations from (1.33|) and (1.36|) we see that (^^) and (2.11) can be rewritten 


Qf{x) 

Qj (x) 

JV/ 






(2.13) 

(2.14) 


It thus follows that the quantities in (2.2) can then be expressed simply in terms of the functions 
and K:^ introduced in ( 1.32 ) and ( 1.35 ). 


Proposition 3. The (fci, k 2 ) point parity respecting correlation for the PDF (1^ ) is given by (1.24) with 


Kto{x,x') 

Kte{x,y) 

Keoiy,x) 

Kee{y,y') 


= —e 


— A(x — x ')/2 


dx 


-^Ai.-y)l2^^>y ^ ^A./2 r 

^0 


= K^{y,y') 


(2.15) 


(we have appended the superscripts L on the LHS as notation for the Laguerre case Similarly, 

the {ki,k 2 ) point parity respecting correlation for the PDF (1.16) — the Jacobi case to he denoted by 
appending a superscript J — is given by ( 1.24 ) with 

Kooix,x') = -( 1 ^:^) ' J ^{^-u)^^^K;({x,u)du^ 

Ki{x,y) = -(l-x)^/2A{(i_x)i-^/2if„^(x,y)} 

Keo{y,x) = -(yT^) ^ Xx>y + {l-y){l-xr^^^ j ^{l-v)^^^K(,{v,y)dv 


KL{y,y') = {l-y)Ki{ij,y'). 


(2.16) 


We see from ( ^.15| ) that Kf(,(^x,y) coincides with K(f{x,y) which we know from ( 1.34 ) determines 
the fc-point correlation for the LUE. This property of the even-even correlations is equivalent to the 
statement (1.10). Similarly, the final formula in ( 2.16| ) implies that the /c-point correlation for the even- 
even correlations in the Jacobi case coincides with the /c-point correlation for the JUE with parameter 
value (a, 6) i-^ (a,0) (recall (1.37)). This result is equivalent to the statement (1.18). 
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3 Correlations for decimated orthogonal ensembles with a pa¬ 
rameter 

3.1 Quaternion determinant formulas 

To calculate the parity respecting correlations for (O) and (|1.16|) , one notes that analogous to (iH) they 
have the general structure 


W_Wo{X2j-l)We{X2j) {x j - Xk) det[K{x2j-l j X2k]j ,k=l,... 

j—1 l<j<k<2n 


(3.1) 


Explicitly, we choose K{x,y) as in (2.4), Wo{x) and We{x) as in (2.7), and make use of (2.6). The 
significance of this is that for general PDFs of the form (0, a result of Rains Q gives that the parity 


respecting correlations have the form (1.27), and further specifies the elements of the matrix therein. 
Before stating the latter, we remark that the qdet operation in ( 1.27| ) is well defined on matrices A with 
the self dual property 

0 -1 


A = ZrjA^Z. 


2n 1 


Z2n Ir 


1 0 


(3.2) 


This is equivalent to requiring that AZ 2 n be antisymmetric. On the latter class of matrices the Pfaffian 
operation is well defined, and we have in fact that 0 


qdet A = Pf{AZ 2 n), 


(3.3) 


which for our present purposes can be taken as the dehnition of qdet (in fact the results of are written 
in terms of Pfaffians). 

Let us introduce arbitrary polynomials Rj{x) of degree j, and let us follow and introduce the 
notation 

/ OO 

w{y)K{x,y)f{y)dy (3.4) 

-OO 

where we have set 

Wo{x) = We(x) = w(x) (3-5) 


(this can always be accomplished by changing the definition of K{x,y)), as well as the 2n x 2n antisym¬ 
metric matrix with elements 

Mjk= [ dxw{x) [ dyw{y)(Rj{x)Rk{y) - Rk{x)Rj{y))K{x,y), (3.6) 


where the Rj{x) are the arbitrary jth degree polynomials introduced below (1.27). Then according to 
iQ the parity respecting correlations for the PDF ( |3.1| ) are given by (1.27) with the elements of the 
matrices therein specified by 


foo{x,x') 

foe{x,y) 

feo{y,x) 

fee{y,y') 


YTk=o w{x)Rj{x)M^*{K ■ Rk){x') - YTk=o w[x)Rj{x)M-*w{x')Rk{x') 

E ■ Rk){x') - E• R,){x)M-^^w{x')Rk{x') 

EylEo w{x)Rj{x)M-*w{y)Rk{y) - E^IEo w{x)Rj{x)M-*{iA ■ Rk){y) 
Ej.1Eo(«^ • Rj){x)M-*w{y)Rk{y) -n{x, y) - EI1Eo(«^ ' ■ Rk){y) 


-k{x, y) + Eyfc=o(«‘ ■ Rj)iy)Mjk{K ■ Rk){x) - Eyfc=o(K‘ ’ Rjiy))Mj,!w{x)Rkix) 


2n-l , 


Eifco w{y)Rj{y)M.f^{K ■ Rk){x) 


- Eil=o w{y)Rj {y)M^,^w{x)Rk (x) 


jk'’ 


ErEo(^‘ ■ R,){y)M-My')Rkiy') - ■ RkW) 

[ EZ~=lw{y)R,{y)MZw{y')Rk{y') - E ' ^fc)(2/') J 


(3.7) 
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Let us now consider the parity blind correlations. For this purpose we write the PDFs (^) and 
(|1.16) in the form 


2n 


{xj - Xk)Pi[e(xj,Xk)]j,k=i,..., 2 n 

j—1 l<j<k<2n 

where w{x) is given by (^) with the substitution (^), and 


(3.8) 


e{x,y) = 


gA\x v\/'^sgii(^x — y) 

'1 — X\-^ sgn(a:-y)/2 

y 




sgn(x - y), 


Laguerre case 
Jacobi case 


(3.9) 


In (3.9) the notation sgn(x) denotes the sign of x and is thus equal to 1 for a; > 0, to 0 for a: = 0 and to 


— 1 for a; < 0. The equality of (|3.8D with (0) and ( |1.16| ) follows from the identity 


Pf 


l-V/(Xfc)/ 


sgn(a:j - Xk) 






i=i 


f{xQ{2j}) 


(3.10) 


where 


XQ( 2 j-l) > XQ( 2 j), Q{2j) > Q{2j - 1), (j = 1,... ,n) 

and e{Q) denotes the signature of the permutation Q. The identity (|3.10|) can be seen to follow from the 


definition 


pf^ = E 


P{2l)>P{2l-l) 




P{21-1), 


1=1 


valid for any 2n x 2n antisymmetric matrix A, where the * denotes that only permutations which give 
rise to a unique product of the Ojfc’s are to be included. There are (2n — 1)!! such permutations, one of 
which for the Pfafhan (3.10) contributes the term on the RHS of ( 3.10| ). All other permutations give a 
contribution which cancels in pairs, so the Pfafhan ( 3.10 ) is in fact equal to this single term. 

The quaternion determinant formula (1.31) for PDFs of the form ( [3.8] ) has been given by Frahm and 
Pichard |^. In the matrix elements of the 2x2 matrix / are given in terms of skew orthogonal 
polynomials. In keeping with (3.7), we prefer to follow and state the form of the matrix elements 
which involves arbitrary polynomials of degree fc, Rk{x), as well as 


/ OO 

w{y)t{x,y)Rk{y) dy 

-OO 


(c.f. (3.4)) and the inverse of the antisymmetric matrix 

dxw{x) / dyw{y)[Rj{x)Rk{y) - Rk{x)Rj 


/ 

J —( 


ix,y) 


(3.11) 


(3.12) 


Comparing the dehnition (|3.9| ) of t{x,y) with the dehnition (|2.4| ) of K{x,y) we see that the integral in 
( [3.12 ) is unchanged if e{x,y) is replaced by K{x,y), and is thus equal to Mjk as dehned by (3.6). With 
this understood, reading off from |M we have 


f{x,y) = 


Eyfc=o w{x)Rj{x)M-^^ (e • Rk)(,y) - T,jjc=o w{x)Rj{x)M^^^ w(y)Rk{y) 

-<^{x, y) + Ej.fe=o(£ • ■ Rk)iy) - Eifc=o(e • R3){x)M-,yw{y)Rk{y) 


(3.13) 
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3.2 Skew orthogonal polynomials 

Our ability to obtain closed form expressions for (3.7) in the sense of ( ^.15 ) and (2.16) relies on the 
construction of the polynomials {i?j(x)} so that they exhibit the skew orthogonality property 


{R2j,R2k)^ — {R2j+l^ Rlk+l)^ — 0, {R2j,R2k+l) 

where (, )'^ denotes the skew inner product 

ir dy dx - g{y)f{x)), 


= TiSi 




(3.14) 


(/,<?r = 


Laguerre 

' {f{y)g{x)-g{y)f{x)), Jacobi 

(3.15) 

In the case A = 0 the skew orthogonal polynomials in both the Laguerre and Jacobi cases are known in 
terms of an explicit series of Laguerre and Jacobi orthogonal polynomials respectively |Q, as well as a 
more compact form involving derivatives and integrals of these bases Q . Analogous explicit forms of the 
skew orthogonal polynomials are available in the limit A —> —oo, when the PDFs ( |1.2D and (1.16) tend 
to certain PDFs of the form ( 1.15 ) corresponding to a symplectic symmetry. Regarding this latter point 
note that integration by parts shows 


lim 

A—^ —OO 


A\2n 


ii)' 


dxi 


'Xl 


■dx2n ]^a;(x20f 


-(l+A)x2i-i/2 -{l-A)x2i/2 


n 


l=\ 


l<j<k<2n 


fXL 


dX 2 dXi ■ ■ ■ dX 2 n ai{x 2 l)e {X 2 j - X 2 kY, 


(3.16) 


1=1 


1<j<fc<n 


where the ai are arbitrary, Al is the integration region and Xl the integration region X 2 > X 4 > 

• • • > X 2 n > 0, and 

j 2 Tl P 

^lim (^—'j / dxi • • •dx2n J|ai(x2/)(1 - X2z-i)^““^“^^/^(l - {xj - Xk) 


fXj 

P n 

/ (iX2(iX4 • • •£iX2„ ]^az(x2z)(l - X2i)“’''’^ {X2j-X2kp, 

1=1 l<j<k<n 


l<j<k<2n 


(3.17) 


where the ai are arbitrary, and Xj, Xj are the integration regions 

1 > Xl > X2 > • • • > X2n > — 1, 1 > X 2 > X 4 > • • • > X2n > — 1 

respectively. Thus in the limit A ^ — oo (0 reduces to the LSE with parameter a = 0, while ( |1.16[) 
reduces to the JSE with parameters a i—> a + 1 and 6 = 0. We see from (3.15) that in this limit the skew 
inner product takes the form 


lim , I 

A — . —oo V 2 / 




J^dye ^/^(/( 2 /)^(e ^/^ff(2/)) - (/^ £/)), Laguerre 

III dy (1 - 2/)(“+^^/^(/(?/)^((l - 2/)(“+^^/^ff(?/)) - (/ ^ 5)) , Jacobi 


In keeping with the known results in the case A = 0 and A —> — oo, we find that the skew inner 
product admits skew orthogonal polynomials with compact expressions in terms of classical Laguerre and 
Jacobi polynomials. Their derivation relies on a special integration formula for the Laguerre polynomial 
Lfe(x) and Jacobi polynomial pjl’^\x), and the latter in turn rely on differentiation formulas for the 
same polynomials. 
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Lemma 4. We have 


p-i 


dt 


Lp{t) — ^ ^ Li{t), 


1=0 


p-i 




1=0 


Proof. The formula (3.18) is an immediate consequence of the well known formula 

■(L“(t)-L“+i(t)) =L“(t) 


£ 
dt' 


in the case a = 0. To derive ( 3.1£ ), we make use of the general formula 

= ^(P +a+ b+ 

in the case 6 = 0, then the general formula 

(n + i(a + 6) + 1)(1 - t)F£+^’^)(t) = (n + a+ - (n + l)P£/£t) 

again in the case 6 = 0, to deduce that 


The formula ( 3.19 ) now follows from repeated use of the general formula 

\ ((2n + a + b)P£’’-^Ht) -(n + 
n + a + b \ 


in the case 6 = 1. 


(3.18) 

(3.19) 


□ 


Proposition 5. We have 


pOO 

g(i-^)t/2 / e“*^^“^)^/^Lfe(a;) dx = '^ cjl^Lpit) 

t ^—n 


where 


and 


where 


^kk 


1- A’ 


p=0 


T , 4 / 1 — y4\ P+l“^ 

cL = {-l7-\, -, (p = 0,...,fc-l), 


f- ) 

{i-A7\i + aJ 

V t ^ _r\ 


p—0 


Cfcfe — 


fc + (a — A + l)/2’ 


Cfcp = (2p + 1 + a)ApBk, (p = 0,..., /c - 1), 


(3.20) 

(3.21) 

(3.22) 

(3.23) 


_ r(p + (g — A + l)/2)) _ T{k + [a + A + l)/2)) 

r(p + (a + ^ + 3)/2)) ’ - r(fc + (a - ^ + 3)/2)) ’ 


(3.24) 


Proof. Consider first ( 3.20| ). Multiplying both sides by e and differentiating, making use of 

(^■ISD , gives 

(1_^) fe k p-i 


—Lkit) = —- 




p=0 


p—O l—O 
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Equating coefficients of Lk{t) gives as stated in ( 3.21 ). Equating coefficients of Lp{t) {p < k) gives 


0 = c^j {P<k). 


j=p+i 


Replacing p by p — 1 and subtracting shows 


Sk,p = (P ^ 


(3.25) 


(3.26) 


Solving this for p = k — l,k — 2,... ,0 in order, making use of the known value of c^p, completes the 


derivation of ( 3.21 ). 

Consider now ( 3.22| ). Multiplying both sides by (1 — and differentiating, making use of 

(^■19| ), gives 

k 


-P, 


(a,0) 




P—0 


k p—1 


+ E E + 1 + a)(-l)^-i-^i^(“’°) (t). 

p=o j=0 

Equating coefficients of (t) gives 

k 

—^p.fc = — (p + ~ ^+ E ^ 

i=p+i 

and replacing p by p — 1, multiplying by (2p + 1 + a)/(2p — 1 + o), then adding to the original we obtain 

1, . 1 (P“t”(^ — ^ — 1)/^) (^P “t“ 1 P rt) 


(3.27) 


c ( , I /I , J -r yu, — -r -r u.) j 

-^P.fe- (p+ 2 (« + dl+l)Jcfcp Ckp-l- 


(3.28) 


It follows from ( 3.27 ) with p = k that the value of is as stated in ( 3.23 ). From knowledge of we 
can use (3.28) with p = A: — 1, A: — 2,..., 0 in order to deduce the formula for c^„, p < fc in d^). □ 

An immediate corollary of Proposition ^ combined with the orthogonalities (1.33) and (1.36) is the 
following integration formulas, which have direct use in the determination of the sought skew orthogonal 
polynomials. 


Corollary 1. We have 




0, j >k 
Ckj, j<k 


(3.29) 


and 


J dA(l-A)(“-^-i)/2pj“’°^(A)^ dx(l-a:)(“+^-i)/2pi“’°^(a;) 


0 , 


2“+i 


j > k 
-cij, j<k 


(3.30) 


2j + + 1 

The result of Corollary ^ implies that in both the Laguerre and Jacobi cases, we have identified a 
family of polynomials {pj{x)} (pj{x) = Lj{x) in the Laguerre case, and Pj{x) = Pj‘^’^\x) in the Jacobi 
case) such that 

ajbk, j < k 


{Pj,Pk) = < 0, 


j = k, 


(3.31) 


-Qjbk, j > k 
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for certain Oj, bk- We can use this special structure to construct the corresponding skew orthogonal 
polynomials as series in {pj(x)}, 


Ri{x) = '^aijPjix), ail = 1- 

3=0 


(3.32) 


This is equivalent to finding a lower triangular matrix T = [ajk]j,k=o,i,..., 2 n-i with I’s down the diagonal 
such that 

0 ro 
-To 0 

0 ri 
—ri 0 


T[ajk\T^ = 


(3.33) 


0 Tn-l 

-Tn-l 0 

where on the RHS all elements not explicitly shown are zero. An explicit solution to this problem is given 
by the following result. 

Proposition 6. Let Pj{x) (j = 0,1,...,) denote a polynomial of degree j, and suppose the value of the 
skew produet {pj,pk) factorizes as specified by (3.31). Then for I even, with 

n ll2—l ■p-rZ/2 T 


n 112 — 1 7 

_ _ _ . _ _ li=3 ^2^+1 11^—j-i-i 2ii ■ ^ I /o 1 /Q q/i\ 

ai2j+i — „i/2-i rr^/2 h ’ ~ rr^/^ ■ ’ ■? — 1 (3.34) 


rl/2 


n lfZ—1 -p-rt/Z 7 

and for I odd, with 


aii-i = — 


h 

bi-i’ 


0.211 n/=j+i b2ij.-i 

aij = 0, j < ; - 2 


(3.35) 


the polynomials ( 3.3^ exhibit the skew orthogonality property ( 3.1^ ). The corresponding normalization 
is given by 

X{l-l)/2 = Ol-lbl. 


(3.36) 


Proof. Suppose first that I is even, and consider 

{pj,Ri)^, j<l. 


Since pj can be written in terms of {Rk}k=o,...,h if follows from (3.14) that 

{p„Ri)^ = 0. 

But on the other hand, it follows from (3.32|) and (^.31) that 


i-i 


{pj,Ri)^ = aj ^ aif_tbfj.-bj'^aif^afj,. 


(3.37) 


(3.38) 


Equating ( ^.37 ) and ( 3.38 ), and calling the resulting equation Cj we see that forming 


-Cj H-Cj—1 

Oj Oj-i 


gives the equation 


h 1-2 i-i 

Oj—\ ^^ Oj ^^ 

-2^ H-2^ = 0, 


a?-! 


(3.39) 


fi—O 


^—0 
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valid for j = 1,2,.... Substituting j = 1, 2,... in order in ( 3.39 ), we deduce that 


OLl2k+l = 

Oil2k+2 = 


0,2k 

a2k+l 

b2k+l 

b2k+2 


-Otl2k 


Oil2k+l 


(3.40) 


valid for fc = 0,1,..., //2 — 1. Recalling the normalization an = 1 we see that the recurrences reproduce 

(Hi. 

Consider now the case I odd. Then 


and it follows from this and (3.38) with j = I — 1 that 


1-2 


r'(/-i )/2 


= {oi-iBi — Bi-i ^ 




^—0 


Making use of (3.39|), which remains valid for j < I — 1, we thus have 


1-3 


^(/-l)/2 . bl-2 

= ai i-ibi-i - - - 2_^ ai^a^ 


bi- 


1-2 


Ol-l 


Ol-2 


Ol-l 


-E 




(3.41) 


(3.42) 


^—0 ^—0 

The fact that ( |3.39| ) remains valid for j < I — 1 means that the equations (|3.40| ) are again valid, this time 
for 2fc + 1 < Z — 2 in the first equation and 2fc + 2 < Z — 3 in the second equation. In particular, since we 
are assuming Z is odd, it follows from the first equation in (|3.40) that 




TO = 1, 3,..., Z — 2 


^=0 


and thus ( 3.41 ) and ( 3.42| ) simplify to ( ^3.36 ) and 

r{l-l)/2 


, bi-2 

= Oili-iOi-i - aii-sai-s 


(3.43) 


Ol-l ai-2 

respectively. Substituting ( 3.3(]| ) in ( 3.43| ), we see that choosing q ;;;-3 = 0 implies the value of on-i given 
). The final equation in 


m (|3.at)| J. me tmal equation m ( |3.at| ), aij = 0, j < Z — 2, follows from having chosen an-^ = 0 in 
and the recurrences ( 3.40| ). □ 

Examination of the above proof shows that for I odd the value of ai ;_i is in fact completely arbitrary. 
This is because the skew orthogonal polynomials as specified by ( 3.14 ) are not unique. For a given family 
of polynomials {Rjix)}j=o,i,... satisfying (3.14), the family with 

R 2 j+lix) 1 -^ R 2 j+l{x) +^jR 2 jix), 

arbitrary, also satisfy O- This non-uniqueness underlies the arbitrariness of au-i; the choice made 
in ( 3.35 ) leads to the simplest result in that with this choice we then have aij = 0 for all j < Z — 2. 

Inserting the explicit value of ajbj. in ( |3.31 ) from Corollary |^, we get from Proposition ^ the following 
explicit formulas for the skew orthogonal polynomials in the Laguerre and Jacobi cases. 


Corollary 2. The polynomials 


Rflllix) = L 2 l-l{x)- 


A + 1 
A-1 




1-1 


(2^) = (^) - 4^ (a 


1=0 


(3.44) 


1=0 
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are skew orthogonal with respect to the skew inner product ( 3.11 ) in the Laguerre case, while the polyno¬ 
mials 


2 /-n J 2 i-i \ ) 2l + {a-A-l)/2 ^'-2 ^ '' 


2^_2 


3=0 


21 + {a- A + \)/2 


(3.45) 


are skew orthogonal with respect to the skew inner product (3. H; in the Jacobi case. The corresponding 
normalizations are 


JL) _ 

I 1 — 


(1-A)2’ 


AJ) _ 

'i — 


2“+i 


{2l + {a- A + i)/2){2l + (a - A + l)/2) ■ 


(3.46) 


Even though the matrix elements (3.7) for the correlations (1.27) explicitly depend on {i?j(a:)}, we 
will not directly make use of the formulas ( 3.44 ) and ( 3.45 ) in our subsequent simplification of (3.7). 
Rather we will make use of these formulas to evaluate the indefinite integral 


w{x)k{x, y)Rk{x) dx, 


(3.47) 


which will then be used in With Rk{x) replaced by Pk{x), this integral is given by ( |3.20D in the 

Laguerre case, and ( [1.22| ) in the Jacobi case. Using the notation of the RHS of (3.31), and introducing 
the additional symbol dk, these results can be combined into the single formula 


1 


w{y) 


w{x)K{x,y)pk{x) dx = ^Pk{y) + X! -^Poiy) 


k-1 


dj bk 


j^o 


where 


7(a;) = 


hj := 


w{x), Laguerre 

(1 —a;)r(;(a:), Jacobi 

We can use ( 3.48|) together with the result of Proposition || to evaluate ( 3.47|). 


w{x)w{x){pj{x))‘^ dx. 


(3.48) 


Proposition 7. Let {Rj(x)} be given by (3.3i), with the aij therein specified by Proposition^. Further¬ 
more, assume the integral evaluation ( 3.46 ). Then 


where for k even 


while for k odd 


'Ky) 


Ukl 


Ukl = 


1=0 


w{x)k{x, y)Rk {x) dx = E ukipi {y) 

otkidi aittki+ih+i 


ctkiai 

hi 


hi 


ak 

Vk' 

dkOikk-i , cik-ibk 


I odd 
I even 

I = k 

l = k-l 
otherwise 


(3.49) 


(3.50) 


(3.51) 
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Proof. Substituting ( 3.32 ) in ( 3.47 ) and making use of ( 3.48 ) gives 




/ w{x)K{x,y)Rk{x) dx = ^ + ^akj ^^pAv)- 

V j=o J 


j=0 tj.=0 


=n ^ 


The coefficient of piiy) in the above expression is 


(Xkiai , ai 


j=i+i 


(3.52) 


But for k even, the second formula in (3.40) shows that we get cancellation in pairs in the above summa¬ 
tion, and (3.5C) results. For k odd, we see from ( [1.35 ) that the summation in (3.52) vanishes for I < k— 1, 
and that so too does the hrst term. The terms which remain give ( [l.51 ). □ 

We now substitute the particular values of the quantities d;,6j,afe implied by (3.21) and (3.22), 
together with the explicit formulas for aki implied by (3.44) and ( 3.45| ), and the normalizations (1.33) 
and ( 1.36 ), in ( 3.50 ) and ( 3.51 ). This shows that in both the Laguerre and Jacobi cases the coefficients 
Uki, up to a sign, are independent of 1. Furthermore, in the case k even we can identify the resulting 
series as a linear combination of 


1 

— / w{x)pk{x)dx and Pfe(a;), 

[y) Jv 


w{y) 

or alternatively as a linear combination of 

1 f°° 

/ w{x)pk+i{x)dx and pk+i{x). 

W[y) Jy 

Corollary 3. The polynomials ( 3.44U have the properties that 

j (t) dt = ( - L^k+Ay) + T2fe 

and 

POO O / ^ ^ 

Jy i \ 

1 


(3.53) 


1-A 

1 

1-^ 


1=0 

L2k{y) 


1=1 
ey/2 

2(1-kl) 


POO 

/ e~'^/'^L2k{s)ds 

J V 


pv/^ r°° 

L2k+iiy)- 2{1-A) j e "/^L2fe-Hi(s)ds (3.54) 


while the polynomials \3.4^ have the properties that 


(a,0) 


(3.55) 


and 


(i_y)-4+i)/2 ^'\[i\t)dt 

J y ^ ^ y'' 


2k 


2k + (a-A + l)/2^ ' 
' " 1=0 


1 


2k + {a-A + l)/2 
1 

2k + {a- A + l)/2 


(AtAy) + fo - i)<-«24”’«)<<*) 


(- / 0 - A ( 3 . 56 ) 
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Proof. The series expansions follow immediately upon making the stated substitutions. To obtain the 
integral formulas, we substitute for the summations in (3.54) and ( ^.56 ) according to their value implied 
by ( ^.2(]| ) and ( ^.22 ) respectively. □ 

We remark that the series of Laguerre and Jacobi polynomials in (|3.53)-(3.56) can each, according 


to the results ( ^.44| ) and ( ^.45 ), be identified with Rj{y) 


A —>• —c 


for j = 2fc + 1 or j = 2k as appropriate. 


3.3 Summation formulas — the even-even block 


The skew orthogonality property ( 3.14 ) implies the matrix Mjk as specified by (^) is equal to —1 times 

], {j,k) ^ (21,21 + 1) 01 {21+ 1,21) 


the RHS of (HH). Thus 


M-k = { -rr\ k) = (21, 21 + 1) (1 = 11,...,n-l) 
{j,k) = (21 + 1,21) (l = Q,...,n-l) 


' i 

rT\ 


(3.57) 


and so the double summations in (3.7) all collapse to single summations. In particular, with the entry in 
row s, column s' of the matrix fab (a, b = eoro) denoted , we have that 


n—1 .. 

fL\y,y') = E - ^hiy')^' 

_ V T 


- To 

j=0 ^ 


2j+l 


(3.58) 


where is defined by ( 1.28| ) and ( 1.29| ) in the Laguerre and Jacobi cases respectively. The latter indefinite 
integrals are precisely those occurring in Corollary |[ Using this result, can be expressed in terms of 
the functions K^n in the Laguerre case and in the Jacobi case. 


K. 


(y, t) dtj 


Proposition 8. In the Laguerre case 

1/7°“ 

feeiy, = 4 (y ^Liy\ t) dt - 

while in the Jacobi case 

feeiy, y') = 7 ((1 - y') [ ^Liy', t) dt-(i-y) j Ki^(y, t) dt). 

^ J y Jy' ' 


(3.59) 


(3.60) 


Proof. Consider the Jacobi case. Substituting for r^^p using ( 3.46| ) and substituting for ^ 2 k+i(y') using 
( p.55 ) we see that 

/ie (2/, y') = ^(1 - 2/)(“+'^/^(i - 

n— 1 

X E(2j + (a - ^ + l)/2)(pE°i(2/')‘i>2, (2/) + 4“’°^ (2/')‘&2,(2/)) - (j/ - y')- 

i=o 

In the product P2‘j+\(y')^2j(y) we substitute for (1 — y)^“'''^^/^$ 2 j(y) using the third equality in (3.56), 
while in the product P 2 'j'-^l(y')^ 2 j(y) we substitute for (1 — using the second equality 

in ( 3.56| ). Straightforward simplification and comparison with ( |1.35D , taking note of (1.36), gives (3.60). 
The Laguerre case is similar. □ 

The single sum form of the other matrix elements of fee are 


n— 1 


/ee (2/,2/') = /ee (2/',2/) = ” E ^ (2/)i?2, + l (y') - +,(y) Ry (V )) . 

Ti \ / 


j=0 
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n—1 


feeiy,y') = -Yl 

J=0 


w{y)w{y') 


(i?2i(y)i?2j+i(y') - i?2i+i(2/)i?2j (2/'))- 


(3.61) 


Noting from (1.28) that in the Laguerre case 




and from (1.2S|) that in the Jacobi case 


= -(1 - 

we see that all the quantities fH, and can be expressed in terms of and thus and K^^- 


Proposition 9. In the Laguerre case 

'A d 


feeiy^y') = /ee (2/'>2/) = (f + fee iV^y'), feHv^y') = ^ ^ (3-62) 


dy' 


dyy 


dy' 


while in the Jacobi case 


feHy^y') = feKyfy) = + -^)feeiy^y')^ 


f!eiy,y) (2^_ ) + a„)(2(1 _„') + 


A 


d 


.2(1-2/) ' 52/A2(1-2/0 ' 

The (0, A:)-point parity aware correlation, or equivalently the fc-point correlation for the even labelled 
coordinates, is according to (1.24) equal to qdet[/ee(2/j, yz)]j,;=i,....fc- By performing elementary row and 


column operations, making sure to conserve the self dual structure (^), we see from Proposition ^ that 
in both the Laguerre and Jacobi cases all terms dependent on the parameter A can be eliminated, leaving 
as the final expression 

d 


P(Q,k){yi,---,yk) = qdet 


gyJeHyj,yi) flKyj^yi) 

' f^^ ’dJi ’ 2 ^' ) ^ f^"^ , yj ) 

dyjdyi dyj 


(3.64) 


J j,l=l,...,k 


The independence on A is required by the identity (1.19), valid for the Laguerre and Jacobi weights 


in (1.20). Moreover, this fc-point correlation must agree with the fc-point correlation for the LSE with 
parameter a = 0 in the Laguerre case, and for the JSE with parameters ai—>0+1, 5 = 0 in the Jacobi 
case (see ( |3.106D and Section |5.2D . 


3.4 Summation formulas — the parity blind case 


The matrix element (3.61) is fundamental with respect to all other matrix elements in ( |3.7D and ( |3.13 ) 

in the 


in that each of the latter can be constructed from ( 3.61 ) by integration. The evaluation 
Laguerre case and (3.63|) in the Jacobi case of (3.61|) then allows us to express all matrix elements in (^?7| 


and (3.13) in terms of and In this subsection we will undertake this program for the matrix 

elements of (3.13). Formulas for the matrix elements of the blocks feo, foe and foo in (^1.17 ), obtained 


using knowledge of the evaluation of the matrix elements for the block /ee given in subsection 3.3 and 


the evaluation of the matrix elements of ( ^.13 ) to be given in this subsection, will be presented in the 
next subsection. 
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Recalling the definition of /ee (y,yO from (3.7), and the definition ( 3.11 ) of (e ■ Rk)^ we see that the 
matrix elements of f(x,y) in (3.13) can be written in terms of f‘^1 according to 

(21! 


f {x,y) = r\y,x) = / e{y,t)k^{x,t)dt (3.65) 

J —OO 

f^^{x,y) = -fll[x,y) (3.66) 

/ OO pOO 

dse{x,s) / dte(y,t)fll[s,t) (3.67) 

-OO j —OO 

These formulas can be made more explicit. For this purpose use will be made of the following formulas. 
Lemma 10. We have _ „ 

poo poo 

- K^uiO,u)du= K^„{y,u)du 
Jv Jo 


and 


-J K^ni-^,u)du= J ^K^^{y,u)du 

Proof. It follows from ( |3.20D with t = 0, A = 0, and the evaluation formula 

LpiO) = 1 

that 

/ e-*/^Lp{t)dt = 2{-lf. 

Jo 

Recalling the dehnition ( 1.32| ) of we thus have 

poo 2n—1 

/ Kt{y,u)du = 2e-yl^ ^ {-IfLp^y). 

Jo 

Also, from ( 3.2(]| ) with A = 0 we have 

poo 

gt /2 / e-"/2Lfc(a;) dx = 2Lk{t) + 4(-l)'= ^(-l)PLp(t). 


(3.68) 

(3.69) 

(3.70) 


p—0 


k-1 


(3.71) 


(3.72) 


p—0 


It follows from this integration formula and ( |3.70D that the LHS of (|3.68| ) reduces to the RHS of ( 3.71 ). 
To establish ( 3.69 ), we note that it follows from ( 3.22 ) with t = 0, A = 0, and the evaluation formula 

p/“-°)(0) = (-!)' 


that 


J {1-x)^^-^'>/^P^‘'’°\x)dx = 


It follows from this and (1.35) that 



(3.73) 

2(a+l)/2 

(3.74) 

k {a l)/2 

2n-l 


(3.75) 




Regarding the LHS of ( |3.69| ), we note from ( 3.22| ) with A = 0 that 

(1 - t)-(“+l)/2 J\l - dj. = 1 


k-1 


+ (a + l)/2 


{Pfc(t) + 2^P, 


1=0 


(3.76) 


This integration formula together with (3.73) shows that the LHS of ( 3.69| ) reduces to the RHS of (3.75). 

□ 
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Proposition 11. In the Laguerre case 


f^{x,y) = 2^2nix,y) - *'>/^K^^{t,y) dt 

A PX POO O A POO 

-- / du—K^^{u,t) + -e^^/‘^ K^^{y,u)du 


f\x,y) = 

while in the Jacobi case 

r{x,y) 


- / K^n{x,t)dt\ 

Jy ^ 

_eA|.-y|/2ggn(^ -y)-{J^" e^^y-*^^^K^^ix,t) 

- l\^^^-*'^/^Ktiy,t)dt} 


dt 


1 If a > ^ — f \ Ail 

-(l-x)ii:2“'„(a;,2/)--|(l-y)—-l| J ^[y^J ^Lit,y)dt 

/ 1 _ / \ Ail r ^ 

^(t^i) / {'-('- 


4(1-y) J_ 
A 


iThrj: 


KL{y^u)du 


r{^,y) = 


4(1-2/) Vl- 
1 / A d 
4 V2(l — x) 9xV V2(l — y) ' dy 




= -(1^) 


' 1 — X \ -^sgn(rc-y)/2 

y 


( fV / l-t\A/2 ^ 

sgn(x - y) - |(1 - x) J ^ (y3^) dC^nixA) 


/ 1 - t \ ^/2 , 'I 

-(i-y)y 1 (y3^) 


(3.77) 


(3.78) 


(3.79) 


dt 


(3.80) 


(3.81) 


(3.82) 


Proof. Consider first the Laguerre case. The formula ( ^.78D follows immediately from ( |3.66 ), upon 
substituting ( 3.5E ) in the second formula of ( 3.62 ). In preparation for deriving ( 3.77 ), we note that the 
last substitution, after computation of the corresponding derivatives where possible, yields 


1 . A s. 2 ( P^ ') 

fee{y,y') = 4(2) jy^, Kt{yA)dt + j^ Kt{y',t)dtj 

A ( f) f) ^ 


Thus there are three distinct terms which must be substituted in (3.65). 

Substituting the A independent term from ( ^.83 ) in ( 3.65| ) and integrating by parts gives the contri¬ 
bution 


1 f A r°° 

_ e^“/ 2 if 2 ^„( 0 , y) + 2Kt{x, y) + ^ y) dt 

-2^ £ y)dt+-^ y) dt). (3.84) 
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The expression obtained by substituting the term proportional to A from ( ^3.83 ) in ( 3.65| ) can be simplihed 
by integrating by parts immediately in one term, while in the other hrst making use of the identity 


and then integrating by parts. Doing this allows the contribution to be written 

1 


(3.85) 


^Ax/2 


A 

dy' 


1 f) 

K 2 u{t, y)dt- - j y) dt 

A f) 

+ j J du—K^^{u,t)- 


(3.86) 


Simplifying the contribution to (3.65) from the term proportional to from (3.83) in the same way 
gives for the final term 


A r r°° r°° r> 

_|_gAx/ 2 y K^jO,u)du- du—K^Ju,t) 

poo poo 

K2niy,u)du- 


(3.87) 


Adding together ( 3.84 )-( 3.87 ) and simplifying using ( |3.68 ) and the formula 

/•OO Q 

KLiO,u) = J —K^„{y,u)du, 


(3.88) 


which follows from (|3.68D by differentiation, we obtain ( 3.77 ). To derive ( 3.79 ), we first note from (^ 
that (3.67) can be rewritten 


/ OO 

dse{x,s)f^{s,y), 

-OO 


and then substitute (3.77) in this formula. Simplification along the same lines as that detailed above 
gives ( pJ^ . 

Next we turn our attention to the Jacobi case. Substituting ( [3.6C ) in (3.63), then substituting the 
result in ( 3.66 ) gives ( t3.8l| ). The former substitution, with the derivatives computed where possible, 
yields 


7 L (i-(‘-!''>|^)"'''i(»'.*) 


A ( 1 


1-y' 


dt 

(3.89) 


((1 - - i) - (li - k')^ - i)} i<L{y. y')- 

This has an analogous structure to ( ^.83 ), consisting of three distinct terms which must be substituted 
in (|t^) . 

Substituting the A dependent term from (3.8£) in (3.65|) and integrating by parts gives the contribu- 
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tion 


A / 1 — t\ Asgn(x-t)/2 


l-t\^/2 


1 — X 


K2nit,y)dt 


KL{t,y)dt 


7 (77)' 

)i:m) 




r-l / X — t \ Asgn(x-t)/2 


KL{t,y)dt. 


(3.90) 


For the term proportional to A in (3.89), as well as using direct integration by parts to simplify its 
contribution to (^.65) in one of the terms, we make use of the identity 


_: 7 sgn(a; - t)-^(l - t)^^Snix-t)/2 ^ _ ^yAsgnix-t)/2 

2 1 i di 


(3.91) 


(c.f. (3.85)) and then integrate by parts in the other. This contribution is then found to equal 

-liThT'^L ('“''■ 


1 / 1 — f \ Asgn{x—t)/2 / (9 \ T 


+ 


/-I 

A 1 




(3.92) 


Use of (|3.9l| ) and integration by parts shows that the contribution to ( [3.65 ) from the term proportional 
to A^ in (|3.89|) can be written 


+ - 


1 / 2 n^/2 

f KL{-h 

— y\l — xJ 

Jy 

1 (l-t^ 

24sgn(a: — t) /2 / 

H 

1 

_I 

I 

' ( 


) du 


(/ (l-(l-i)-^)dCL(i,u)duyt 


A/ 2 \A/2 A yl-t\Assn{x-t)/2 

(l^) 8 L(t^) (3.93) 


Adding together ( 3.90D , ( |3.92 ) and ( 3.93 ), and making use of ( [3.69| ) gives ( 3.80 ), but with the term 


-{l-x)K:^^{x,y) 


(3.94) 


replaced by 

/ 2 \ a/2 , ^ / B \ ^ — t \ a/2 

In fact ( |3.94 ) and ( 3.95 ) are equal, as is deduced by integration by parts in the latter. □ 
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In the orthogonal symmetry limit, A = 0, the results of Proposition 11 can be written 
1 \ d 

7 ;K 2 n{x,y)---^ K^^{t,y)dt, Laguerre 

, 5 N r z 

-{l-x)K^^{x,y)--{^{l-y)—-lj J ^Ki^{t,y)dt, Jacobi 



^f^{x,y) = -^sgn{x-y)-J f^{x,t) 


dt. 


(3.97) 

(3.98) 


The form of /^^(x,j/) given here is different to the form known from earlier literature (see Section 5.1). 

It is also of interest to specialize Proposition ^ in the symplectic symmetry limit, A —oo. Inte¬ 
gration by parts shows 


ri^,y) 






Laguerre 


1 19 

2^2;(2;,2/) + 2^y KL{i,x)dt, 

i(l-a:)i4:/„(a:,y)-b^|—^((l-x)^-l)y' K^^{t,x)dt, Jacobi 

- J (t, y) dt, Laguerre 

r P\t,y) 


(3.99) 


, , ^ ■ dt, Jacobi 

1 - a: (1 - t) 


A 

dy 


r\x,y), 

d 


Laguerre 


(1 - y)^(l - y)/^^(a:,y), Jacobi 
dy 


(3.100) 


(3.101) 


Equivalently, these formulas can be written 

l( f K^n iy,t)dt- f {x, t) dt\, 

1/1 " I t 

2\T^ J ^^riiyA)dt-J Ki„{x,t)dty Jacobi 




Laguerre 


A^ 


r\^,y) 


r\x,y) 


^ r2i/ ^ 

{y,x), Laguerre 

{l-x)-^{l-x)^f^{y,x), Jacobi 


d ^2 


f^ix,y). 


Laguerre 


dydx 8 ^ ^ 

-(1 - ai)(l - - x){l - y)^f^{y,x), Jacobi 


(3.103) 


(3.104) 


Defining 


r\x,y) = 


1™ 


A^- 


Laguerre 


A^ 


(3.105) 


(1 —x)(l —y) lim -—f^{x,y), Jacobi 




it follows from ( 3.102 )-( 3.104 ), the first equality in ( 3.65| ) and the quaternion determinant formula ( 1.31 ), 
that the A:-point correlation is identical in structure to ( |1.64| ), but with replaced by Furthermore, 
comparing ( 3.102 ) and ( ^.105 ) with ( 3.59 ), ( ^.6C ) we see that 


f^^{x,y)=2f^,^{x,y). 


(3.106) 
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This is consistent with the fact that both the even-even correlation ( 3.64 ) and the A —oo limit of the 
parity blind correlations coincide with the correlations for the corresponding symplectic ensemble (the 
factor of 2 in ( ^.101^ ) is due to the double degeneracy inherent in the A —> —oo limit). 

3.5 Summation formulas — the odd-even and even-odd blocks 


The blocks feo and foe in (3.7) are duals in the sense of (|3.2D, and thus 

feo{y,x) = 


foe(.x,y) -ff^{x,y) 
-foeix,y) foeix,y) 


Consequently, it suffices to consider one of these blocks, foe say. Now, analogous to the formulas ( |3.65 )- 
( p.67D expressing the elements of the matrix f{x,y) in terms of fff, we can express the elements of 
foe{x,y) in terms of fff. Thus we see from 0 that 

fllix.y) = feeix,y) 

pOC) 

foeix,y) = - K{t,y)foe{x,t)dt = -/“(y,x) 

'^y 

folix,y) = f KixA)feHt,y)dt 
J — oo 

/ X poo 

dtK{x,t) / dsK{s,y)fll{t,s). 

-OO J y 

The last two formulas can be rewritten to read 


(3.107) 


y) = -f^{x, y) + fll{x, y) 

poo 

fll{x,y) = -K{x,y)- / K{s,y)f^{x,s)ds + fl^{x,y). 

J y 


(3.108) 


The only quantity in these formulas which is not known explicitly from our study of the even-even blocks 


and the parity blind case is the integral in ( 3.1081) . Its simplihed form is readily computed. 
Proposition 12. In the Laguerre case 

poo 1 /■“ ^ r! 

J K{s,y)f'^{x,s)ds =2] - J dt y du—K^^{u,t) 

1 

+ _e^-/2y Kt(0,t)dt (3.109) 


while in the Jacobi case 

roo 


/ 


1 / 1 — f \ 

{s,y)f^{x,s)ds = -{l-y) J ^ 2 n{i,y)dt 

1 7^ /I—7^ / d\ / 2 \ ^/2 7^ 

2 y j (t^J j ^Li-^A)dt{3A10) 


Proof. In the Laguerre case P^{x, s) is given by ( 3.77 ), while in the Jacobi case it is given by ( 3.80 ). 
Substituting these formulas as appropriate, and simplifying according to the strategy of the proof of 
Proposition II gives the stated formulas. □ 
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3.6 Summation formulas — the odd-odd block 


We read off from (H) that 

fooi^^x') = 


foo{x,x') = -f^^{x',x) 

/ X px' 

dtK{x,t) / ds k{x , s) f^l{t, s) 

-OO J —CXD 

= f^ix, x') + fl^{x, x') - /oe (x', x) - fll{x, x'). 


(3.111) 


The quantity fH is specified by ( ^3.10^ ) (with and therein having the explicit forms ( ^.77| ), ( |3.80 ) 
and ( ^.62 ), (3.63) respectively); is given by the explicit forms (3.79), (3.82); by ( ^.59 ), (3.6C); 


foe by (3.1081) (with fff given as noted, and the integral by Proposition |l2|). 


4 Scaled form of the correlations 

4.1 Superimposed orthogonal ensembles with a parameter 

Consider first the Laguerre case. To leading order, the support of the spectrum for the LOE is [0,4n]. 
As has been identified in previous studies (see e.g. |I^), there are three distinct scaling regimes in which 
different limiting forms of the correlations are obtained. These are the hard edge, specified by the change 
of scale 

(4.1) 

(4.2) 

(4.3) 


Xj 1 -^ XjlAn, 

the bulk of the spectrum, specified by 

Xj c + ttXj / y/n, 0 < c < An (fixed) 


and the soft edge, specified by 


Xi 


An + 2{2n)^/^Xj. 


In general, under the linear change of scale 


= a{n) + b{n)Xj, 


the correlation functions transform to correlation functions in the new variables {Xj} according to 

PkiXi,.. .,Xk) = {b{n))^pk{xi,.. .,Xk). (4.4) 

The significance of the particular scales @-(0) is that 

^scaied(^i^..., Xk) := lim {b{n)fpk (a{n) + b{n)Xi,a{n) + b{n)Xk) (4.5) 

is well defined. 


For the parameter dependent extension of the LOE (1.4), and the parameter dependent extension of 


the superimposed ensemble LOE U LOE ( |l.8[) , we expect the hard edge, bulk and soft edge scaled limits 
to again all be well defined provided the scale of the parameter A is suitably chosen. The correct choice 
can be anticipated by the requirement that the quantity occurring in the formulas ( 2.15D and 

(tl.77)-(3^) be of order unity in the scaled limit. This is achieved by 


4na, hard edge 

A 1 -^ < y/najir, bulk 

a/2(2n)^/^, soft edge. 


(4.6) 
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where a denotes the scaled parameter. 

Let c, w > 0 be otherwise arbitrary hxed real numbers. Using the asymptotic formulas 


g-./2^a/2^a(^) ^ „/-a/2 (2(n'x) ) + i?i 

n\ 


(4.7) 


n' — 71 + (a + l)/2, i?i — 
where Ja{z) denotes the Bessel function of order a 


a;5/4o(f^“/2 3/4^^ ^ < X < UJ 

x“/2+20(na), 0 < X < cn~^ 


g-./2^a/2^a(^) ^ ^a/2 1 cos(2(7ix)-0^2"7r/4) + (nx)(1)^ , cn-1 < X < w (4.8) 


and 


e-“/2x“/2L“(x) = ti“/2 


(- 1 )^ 


-Ai(t) + 0(e-‘)o(7i-i/3) 


^2“(2n)i/3 

where x = 47i + 2 + 2{2n)^^^t, t G [to, oo), it is straightforward to derive the well known formulas 


(4.9) 


^liard(^^y) _ 


lim 

n^oo 4 n \471 4 ?! 


= XX,Y>Q 


Ja(Al/2)yl/2j'(yl/2) _ Xl/2 J'(Al/2) J^(yl/2) 


2(A - Y) 

2 ^ rl 

Xx,y >07 / Ja{'/xi)Ja{VYi)dt 

4 Jo 


a=0 


^bulk(^^y) _ 


/o 

TT / ttX -kY \ 

hm —=K^{c-\ -^,cH-^ 

\/n \ \/n ^n) 


a=0 


n^oo Un 


sin7r(A —F) 

- ——I cosTr{X-Y)tdt 

7r(A — y) Jq 


(4.10) 


(4.11) 


^soft(^^y) _ 2{2ny/^K^{An + 2{2ny/^X,4n + 2{2nf/^Y) 

n —>-oo 

Ai(A)Ai'(F) - Ai(y)Ai'(A) 


X-Y 


poo 

/ Ai(A + t)Ai(F + t) dt. (4.12) 

JO 


We see from (2.15) that knowledge of (4.10)-(4.12) immediately gives the scaled form of K^^{y,y'), since 
we have 


^scaied(y^y/) _ liui b{n)K^{y,y') = K^‘=^'^^{Y,Y'). (4.13) 

n —*oo 

The scaled form of the remaining quantities in ( 2.15| ) is also obtained by formally replacing AT^ by if scaled ^ 
We will see that the derivation is straightforward in the cases of the hard edge and bulk limits. For the 
soft edge limit the integrations in the formulas for AT^ and are to leading order over the interval 
[Ojdri] whereas (4.9) applies to the interval [47i,oo). To overcome this difficulty the following identity will 
be used. 

Lemma 13. We have 

J\-^-/^Kf:iy,u)du = e-(i-^)“/2(|-A„(w)) dw 

poo 

s^^/2 / e-^^/^K’i;{y,u)du. (4.14) 

J X 


—e 
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Proof. It follows from the integration formula |2^ 


that 

pOO cy n—1 . _ . , 

I e-^“/2if^(y,u)du= (^) Li{y). (4.15) 

But the sum in ( 4.15| ) can, according to ( ^.2C ), be written as an integral involving a single Laguerre 
polynomial, and ( 4.14 ) results. □ 

Proposition 14. Let the term “scaled” refer to any of the hard edge, bulk or soft edge limits as specified 
by (4.1)-(4-t). Then in addition to (4-15) we have 


Alibied(y, X) 

:= hm Kt,{y,x) = -e^ 

n —>-oo 

^scaled y) 

:= Van {b{n)fKff^{x,y) 

n—^oo 

RScaled^X, X') 

:= lim b{n)K^^{x,x') = 


■■u 




(4.17) 


rX' 


(4.18) 

valid for a < 0 in the hard edge and bulk cases, and for all a in the soft edge case. 

Proof. Using the differentiation formula 

lL;it) =-L;flit) 

one can check the well known fact that the asymptotic formulas (|3)-(|;|) remain valid after differen¬ 
tiation. This, together with ( |4.13 ), then implies (4.17). Consider next The explicit form of the 

remainder term in ( [4.7|) implies it does not contribute to the scaled limit ( [4.l[) , rather the sole contribution 
comes from (4.13D, and this implies (4.16) in the hard edge case (note that the lower terminal in (4.16) 


(4.19) 


can be replaced by 0 in this case). For the bulk limit we note from ( |4.8| ) and (4.9) that we have 

fX 

X) = y) dv + R, 

J — oo 


where 


i? = lim lim [ e-^^”/‘^K^{v,c) dv, 

e^On^oo \/n Jq 

and in which we are free to choose e < c. Now it follows from that 

\Kl4{v, c)| < Vnf{v), 0<v<e 


where f{v) is integrable. Thus with A < 0 the remainder term R in ( 4.1£ ) vanishes. 

The remaining case is the soft edge limit. We make use of the formula (4.14) and the asymptotic 
expansion (4.9), which together imply 

Jim 6(n)e““/2 / e-°‘”^‘^K^iv,y) dy 

soft edee 0 


^ g-aU24ga(X-V)/2 


pOO pOO 

/ Ai(u) du - {u, Y) du. 

Jy Jx 
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The form ( 4.16 ) now follows after making use of the easily verihed integration formula 

/ OO 

e'^‘Ai(t) dt = (4.20) 

-OO 

and the integral form of in ( [4.12[ ) to deduce that 

/ OO pOO 

e-““/2A:s°f‘(w, Y) du = e-“'/24e-“'^/2 / Ai(t)e“*/2 dt. 

-OO Jy 

Finally, the result ( 4.18| ) follows from the fact, already noted, that the key asymptotic formulas 
( [4.7[ )- (|4.9| ) remain valid under differentiation, together with the method just used to derive ( 4.16 ). □ 

We now turn our attention to the Jacobi case. Here there is a hard edge scaling limit in the neigh¬ 
bourhood of both X = —1 and a; = 1, as well as a bulk limit. From previous studies (see e.g. [^) we 
know the appropriate scales are 


a: 1 — 


A 

Y?' 


a; 1 -^ — 1 + 


A 


for the hard edge atx = l,x = —1 respectively, and 


X X 

X ^ cos00-sin00 ~ cos(0o H-), 0 < 0o < tt 

n n 


in the bulk. For the scaling of the parameter A we choose 


In^a, hard edge at a; = 1 

{ a, hard edge at a; = —1 

na(l — cos0o)/sin0o, bulk 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


This is suggested by the criterion that the term ((1 — a:)/(l — tend to a non-constant order one 

quantity. 

A rigorous analysis of the scaling limits for the JOE and JSE has recently been undertaken by Duenez 
[p^ . Following the methodology therein, we make use of the asymptotic formulas 

P^“’^^(cos0) = ( 7 rn)“^/^^sin-j ^cos-j cos(n'0-I- 7 )-I-Ei, 

n'= n-I-(a-I- 6 -I-l)/2, 7 =— — (a-I-1/2) 

for 0 < 0 < TT, where 

El = uniformly for c/n < 0 < TT — e (c, e<Cl) 

and 

(sin^) (cos^) P,(“’^)(cos0) 
for 0 < 0 < TT, where 

E 2 — 


+ 0 + 1 ) 


n! 


sin0 


Ja{n'9) -\- E 2 


(4.25) 


9^l'^0(n ^/^), c/n<9<'K — e 

0“+2O(n“), 0 < 0 < c/n, 

again uniformly in 0. The asymptotic form ( [l.2l| ) is dominant in the bulk, while (4.25) is dominant for 
the hard edge at a: = 1. To study the hard edge at a; = —1, we make use of the fact that 

cos9) = P///’°-\cos9), 

and then use ( |l.25 ). The strategy of the proof of Proposition 14 then yields the following forms for the 
matrix elements determining the scaled correlations. 
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Proposition 15. The results of Proposition 0 hold for the scaled limit of the matrix elements ^2. 1 in 
the bulk and for the hard edge at x = —1. For the hard edge at x = 1 we have 


K, 


hard / 



1/2 

^hard| 

\Y'J 




S -O 

,/2 

\Y 

) 

Xa;: 


d 

r yT 

yl/2 

dX 



a/2 

d r 



axi 

. Here 

; n 


yl/2 roo 


IX' 


(4.26) 


valid for parameter values a < 3/2. Here refers to but without the restriction a = 0. 


4.2 Decimated orthogonal ensembles with a parameter 

The scaling limits introduced in the previous subsection in relation to the superimposed orthogonal 
ensembles with a parameter all carry over to the decimated orthogonal ensembles with a parameter. 
Furthermore, the strategies used to derive the scaled limits in Propositions 0 and again suffice to 
derive the scaled limits of the matrix elements determining the correlations for the decimated Laguerre 
and Jacobi orthogonal ensembles with a parameter. In the interest of economy of space we will restrict 
ourselves to presenting the explicit form of the scaling limit for the parity blind correlations only. 

Proposition 16. Let the term “scaled” refer to any of the hard edge, bulk or soft edge limits for the 
parameter dependent Laguerre ensembles as specified by 0-0- Then the scaled form of the matrix 
elements ( 3. 7''t )-( 3.7^) have the explicit form 


1 


2ni —2 

X 


lim b{n)f^{x,y) 

n —*-oo 

fsLediX,Y) := lim (b{n)ff^{x,y) 

n—*-oo 2n\-^n 

= I (f+ iw) (I +{ x” *> - r "‘I 


(4.28) 


/scaled(^.^) = lim f ix,y) 
n —>^oo 


2ni—»-n 


= -e“l^-^l/2sgn(X - Y) 


^ J —OO j — <X> ^ 


(4.29) 


Proof. The only new feature of the form of the expressions in ( ^.77 )-( 3.79| ) relative to those in (01) 
is the need to analyze 

pOO 

/ K^„{y,u)du. 

Jo 

This is done by noting from ( |3.71 ) and (3.72) that 

POO -| poo 

J K2niy,u)du =2] e~‘”-/'^L2n(u)du-e~y/'^L2n{y) 

i [ e““/^L2n('«) dw - e“*^/^L2n(y)- 

^ JQ 


= 1 - 


Apart from this, we follow the strategy used to deduce the scaled limits in Proposition |l4|. 


(4.30) 

□ 
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In the Jacobi case the analogue of ( 4.30 ) is afforded by (|3.74D -( ^.76|) . Using the resulting identity, 
and employing the strategies used to derive Propositions E and 1^ gives the following result. 

Proposition 17. The results of Proposition E hold for the scaled limit of the matrix elements ( \3.8(\ )- 

r the ha 

f^iAX,Y) = 


( 3.8< ) in the bulk and for the hard edge at x = —1. For the hard edge at x = 1 we have 

1 /X\i/2 


1 r) /-oo 


poo 

/ t^°‘-^'>/^K'^^^^{t,Y)dt 

J X 


a 




(4.31) 


m,^iX,Y) = 


1, 

r a 

d X 

/ a 

d X 

4' 

~ 

dxJ 

\W~ 

dYJ 


t-^/^K'^^'^'^(Y,t)dt-X^/^ t-^/^K^‘^’^’^{X,t)dt'^ (4.32) 


fil^{X,Y) = -(-) 

Xl/2 


X \ -asgn(X-Y)/2 


sgn{X - Y) 


valid for parameter values a < 3/2, and where X^ard iq (^.IC), but without the restriction a = 0 

(both (/.5i) and (4-3S) have been multiplied by —1 — an operation which leaves qdetffhs^rd] unchanged 
— so as to formally conserve the relations (3.97) and (3.98)). 


5 Discussion 


5.1 The orthogonal symmetry limit 

In the orthogonal symmetry limit A = 0 the matrix elements determining the parity blind correlations 
are specified by ( 3.96D -( 1^ ). We have already remarked that the value of f^^{x,y) implied by ( |3.96 ) 
differs in structure to that known from previous literature ® 0i- To present the latter, which are 
valid for general a > —1 in the LOE and general a, 6 > —1 in the JOE, we introduce generalizations of 

(B , (|i^) , 


n—1 


Yjx^y) := (a;y)“/2e-(-+^)/2^^Lr(a;)Lr(y) 

7^0 Ka 

,(a-l)/2, 


(5.1) 


■= {A-AA-y)y " [A + AA + v)) ^ ( 5 - 2 ) 

1=0 ^l,a,b 

In (|5.1D, L( denotes the Laguerre polynomial of degree I with orthogonality property 


^ r(n + n + 1) 


pOO 

/ Pe-*L):^{t)Lf{t)dt = h() Sm,n, r( ’ 

Jo r(n+l) 

while in (^), P^°'P denotes the Jacobi polynomial of degree I with orthogonality property 


(1 _ t)-(l + tpp(a,b)^^^p(a,b)^^^ ^ hi,J^,n, 


A 

^n,a,b 


r(a + 1 + n)r(6 + 1 + n)2“+''+i 
r(n + l)r(o + 5 + 1 + n){a + 5 + 1 + 2n) 


(5.3) 
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The case a = 0 of (O) and 6 = 0 of (^) reduces to ( 1.32 ) and ( 1.36| ) respectively. 

In terms of (5.1) and (p^), the results of give that the A:-point distribution for matrix ensembles 
OE„(x“/^e“'^/^) (LOE) and OE„((l + a;)^/^(l — (JOE), n even, as specified by ( 1.11 ) are given 

by ( 1.31 ) with the 2x2 matrix / in (|l.3l|) having as its top left entry 

/ r \ 1/2 1 r°° 

f^x,y) = (-) Ll+_\{y) / sgn{x-u)Llt\{u)u‘^/^e-'^/^ du (5.4) 

in the Laguerre case, and 


fii(x w) - (l-xi^i^V^V {xv)\ ^ (a + 6 + n)(a-l + n) ( 

/ {x,y) - (1 + („ + 5 + 2n-l) ^ 

x(l + 2/)*'/2pt’^')(y) 1'^ sgn(a; - du, (5.5) 

in the Jacobi case. The other entries in the matrix / are related to by 

f^{x,y) = f^y,x), f\x,y) = -^sgn{x-y)-J f^{x,u)du, f^{x,y) = ^f\x,y), (5.6) 


x(a-l )/2 


which are identical to the structure of the formulas obtained in Section 3.4 provided we interchange 
and (recall the first equality in ( ^.65 ), ( 3.97] ), ( ^.9S| )). 

We know from ® 01 that (|5.4|) can be rewritten to read 


where 


F,\y) = 




f^{x,y) = K^^^{x,y) + F^{y)F^[x) 


-y^/^e-y/^Ll+_\{y) 


F^{x) = ( - 4x“/2e-"/2L“_^(^) _ sgnix - du^ 

poo 

= -n sgn(x-t)p/^-^e-‘/2(L“(t)-L^i(t))dt. 

Jo 


The formula (5.5) in the Jacobi case can similarly be rewritten. 


Lemma 18. The formula (5.6) has the alternative form 

f^(x,y) = (1 - x)K;[,^ f,{x,y) + Ff {y)F^(x) 

where 




a + 6 + n 


(5.7) 


(5.8) 


(5.9) 


(5.10) 

(5.11) 


F^{x) = ( - 4(1 - x)^‘^+^F2(i + a;)V2pK6)(a.) 

+ (a—1 + n) J sgn(x — u)pj^° 22 ^^\u)(l — ^ n'jb/2 

= _|_ ^ J sgn(x - u)(^{n + a + b)P^°‘F(u) + (n + a) (u)"^ (1 - ■u)^““^^/^(l + u)’"^‘^~^du. 

(5.12) 
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Proof. It follows from 


that 

1 


(g “t” 5 “t” 1 -t- ^)(g b ‘2 2/) 1 


Ha,b+i 2{b+l + l){a + b+l + 2l) 
Using this, and the Jacobi polynomial identity 


(Z + ^ + 1)(1 + = {l + b+ (x), 


we see that 


(1 + ^) E TJ -P/“’'+'’(a:)P/“’'+'’(y) 

^La 


1^0 l,a,b+l 

n— 1 

■ " 'H,a,b ^ 


;=o 


G P 6 P 1 P 2/ 


+ 6 + 1 + i)P/“'‘*'’(s) + (O + oPi 


(a, 6 + 1 ) 


^n-+a,6^0 + ^- 1 + 2?^^ 

Simplifying the term in the brackets { } using the Jacobi polynomial identity 

(a P 6 P 1 P P (a P l)P["t\x) = {a + b + 1 + 21)pI"'^\x), 

recalling the definition (^.2D , and substituting the result in ( |5.5| ) we deduce ( 5.10| ) with F^ix) given by 
the first equality in ( 5.12 ). The second equality in ( 5.12 ) can be deduced from the first by verifying that 
both expressions agree at a: = 1, and then verifying that the derivative of both expressions agrees. □ 
We now equate the RHSs of the final equation in (^.6[), with the substitutions (0) and (|5.10| ) (and 
n ^ 2n), and (3.97), with the substitution ( 3.96| ). In the Laguerre case this gives 

1/9 d 


2 (^ + -^^2n{x,y) = C2n,o(l)L,oi^)^2u,oiy) 


(5.13) 


where 


2^n-l,a 


<a(^) = W2e-“/2L:+l(a;) = -W2-ie-^/2(nL“(a:)-(nP«)P^i 


(5.14) 


(the parameter a has been kept general in ( 5.14 ) for later convenience). Note from ( 5.14 ) that the RHS 

(5.15) 


of (5.13) is symmetric in x and y, as is the LHS. The substitutions in the Jacobi case give 


~ ~ ~ ~ y)^2nix,y) = cin,a,o(t>2n,afii.x)i’L,a,o{y) 


where 


2hi 


-( 


— l,a,f) 


f a + b + n \ 

( ) 

\ 2,71) 

'y2n + a + bJ 


<a.6(:^) = (l-x)(“+l)/2pi“'J+'’(x) 

= (l-x)(“+i)/2(l + :,)U2-i^-^ 

^ ^ ' 2n + a + b 

<a,biy) = -(1 - + a P b)Pi^’^Hy) + in + a)Pt’J)(y)) . (5.16) 


|((n P b)Pi°2'’^{x) P nPi‘^’^'> {xij 
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As with ( 5.13 ), we see from ( 5.16 ) that the RHS of ( 5.15| ) is symmetric in x and y, as is the LHS. 
The identity (5.13|) is the special case a = 0 of the identity Q 


+ -^)^nAx,y) =C^.a(</'n.a(a:)V'n,a(2/) +0n.a(y)V’n,a(a;))- 


A simple consequence of (5.17) is the integral representation 

y) = (^n,a( 2 ; + t^n,a{y + i) + (t>taiy + ^)V’n,a(a^ + O) dt. 


(5.17) 


(5.18) 


The result (|5.15|) suggests an analogue of (|5.17| ), and (|5.18| ), in the Jacobi case. 


Proposition 19. With . specified by and ci ,(t>i Ax), Ay) by ^5.1^, we have 


and consequently, for 0 < x,y < 1, 




/»1 

4a:yAl;[,a,b(l-2a;,l-2y) = (^(j)i^^^^{l-xu)f}i^^^^{l-yu) + (j)i^^^f,{l-yu)f}:i^^^^{l-xu)j-^ (5.20) 

Proof. We generally follow the strategy used in Q in the Laguerre case. But before doing so we note 
the general identity 

(d - X)A + (1 _ „) 3 ) _ ,, 8 + (1 _ „) 8 

V ox oyJ \ x — y / x — y \ ox oyJ 

To make use of this identity we note from the Christoffel-Darboux formula (see e.g. |3^) that 

KM - (d - .)(i - „)).-./dd . .,d . 

X y 


where 


Thus 


a„ = 


r(7r a b 1) 


20+6 _|_ _|_ ^,)( 2 n + a + b)' 


(1 — x)^/‘^{l — y)^!"^ _ 


x-y 


aq(l-x)^ + d-s)^ 


x)(in-i{y] - qniy)qn-i{fi (5.21) 


where 


qn{x) ■.= {l-x)^/\l+xf/^Pi^^'>\x), 

SO the task is to compute the action of the operator on the RHS of ( ^.21 ). 

For this purpose, we note from suitable Jacobi polynomial formulas that 

{1 - x‘^)q),(x) = {ao + aix)qn{x) + AoPn-iix) 

(1 - x‘^)q'„_^{x) = --foqnix) - (cxq + aix)qn-i{x) 


with 


— of 

0(0 H” 0 .\X = -“7 

2(2n + a + b) 
_ 2(n + a){n + b) 

“ 2n + a + 6 ’ 


(fZfi a b) 

2 ^ 
2n(n + a + b) 
2n + a + b 


(5.22) 
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Equivalently 


(l-x)q;(x) 


A(x) 

B{x) 

qn{x) 

_ - x)qn-i{x) _ 


. -C'(x) 

-A{x) _ 

qn-l{x) 


where 


A{x) = ^^±^, B{x)= 


1 + a: 


1 + X 


C{x) = 


7o 

1 + a; 


[qn{x) q„_i(x)] 


Introducing the matrix formulation 

qn{x)qn-i{y) - qn{y)qn-i{x) = 
a straightforward calculation using ( ^.23| ) shows 

^^((1 - x)-^ + (1 - y)^) {qn{x)qn-l{y) - qn{,y)qn-l{x)^ 
[qn{.x) qn-l{x)] 


0 1 

qn{y) 

-1 0 

_ qn-i{y) _ 


1 

1 



X — y X — y 


qn{y) 

A{x) - A{y) B{x) - B{y) 


qn-i{y) 

1 

B 

1 

B 



1 [g«(x) g„_i(x)] 

7o 

ao — OL\ 

qn{y) 

(l + x)(l + ?/) 

1 

P 

o 

1 

P 

00 

_ qn-iiy) _ 


1 


(l + x)(l + 2/) 


(1 - _ y)“/2(l + y)b/2 


2n + a + b 


n{n + a + + {n + a){n + (y) 

. 1 (&2 _ ^ (2n + a + 6)2)(Pi“-^)(x)Pi“jf (y) + Pi“’'’)(y)Pi“jJ)(x)) 


Substituting ( ^.24 ) in ( 5.21 ) gives ( |5.19| ). 

With ( ^.19| ) established, we make the changes of variables 


X = 1 — 2e 


y = l-2e-\ 


0 < s, t < oo 


which gives 


- 2e-^ 1 - 26"*) 


(5.23) 


(5.24) 


= <a.6 (<a.b(l - 2e-*X„,,(l - 2e-‘) + <,,,(1 - 2e-‘)X,,(l - 2e-«)). (5.25) 

A simple calculation using integration by parts verifies that 
4e-(^+‘)<,,,(l - 2e-^ 1 - 2e-‘) = 

X {Ka,ba - 2e-^-^)X,,(l - 2e-‘-^) + <,,,(1 - 2e-‘-^X„,,(l - 2e-^-^)) dz.(5.26) 

satisfies (5.25). The general solution of ( ^.25| ) is ( 5.26D plus an arbitrary function of s — t. But we see 
from the dehnition of ^nab ^'^d ipn a b both sides of (5.26) vanish in the limits s,t oo, 

s — t hxed, so this arbitrary function must be chosen to be the zero function. Putting e“® = x, e“* = y 
in (5.26) and changing variables e~^ = u gives ( ^.20 ). □ 

The orthogonal symmetry limit of the scaled matrix elements of Propositions and is obtained 
by setting a = 0. We see from the results of Propositions ^ and (with and interchanged for 
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convenience) that the fc-point distributions are then given by (1.31) with the 2x2 matrix / having as its 


top left entry 


fLiUx,Y) = - i A 


rX 


K 


scaled 


(t, Y) dt 


(5.27) 


in the case of the soft edge and bulk as given by ( [l.l2 ) and ( [l.ll|) respectively), and 

, V \ \/2 Ir) 

flLUX,Y) =-[-) K'^-^\X,Y) + -—Y^/^ t-^/^K'^^^\t,Y)dt (5.28) 

in the case of the hard edge with singularity as a; —> 0 (/ilhard jg gjyg^ by ( 4.10| ), without the 

restriction a = 0). The remaining entries in / are related to as in (|5.6| ). 

The formula (|5.28D has been read off from Proposition]^, deduced from the scaled limit of the Jacobi 


case. The hard edge singularity as a; ^ 0 is also contained in the results of Proposition ^ 

for the scaled Laguerre case. In this case we read off from Proposition [T^ the alternative formula 

rX 


fLiUx, Y) = r)U=o I Y)U=odt, 


(5.29) 


which differs in form to ( 5.28 ) with a = 1. Our first discussion point regarding the form of /Jcaied exhibited 
by ( 5.27 ) and ( 13.28 ) is to reconcile the two seemingly different expressions in the case of the hard edge 
singularity as a: ^ 0. For this purpose we make use of the following recurrence. 


Lemma 20. Let F) be given by {4-1^), without the restriction to a = 0. We have 

\ Y / a —>^a+l Z 

where it is understood that X,Y >0. 

Proof. We first make use of the Bessel function identity 

Uj'^{u) = aJa{u) - uJa+l{u) 


to rewrite the hrst equality in (4.10), without the restriction to a = 0, as 

^hard.^ _ XY^J,+ ^{XY^)UYY^) - FVVa+i(yl/^) Ja(XV^) 

^ ^ 2{X-Y) 

The stated formula now follows by replacing a by a + 1, using the Bessel function identity 

MJa+2(u) = 2(a + 1) Ja+l(u) - uJa{u), 

and simple manipulation. 

Using Lemma ^ in ( 5.28|) we obtain 


□ 


fsLUx,Y) 


a^a+1 2 4 


1 r) 1 f) P^ 

^2WJx X'^^’^''it,Y)dt--^MY^/^) t-i/Va+i(ti/2)dt.(5.30) 


In the special case a = 0 the last integral can be evaluated, and this shows the final term cancels with 
the second term. We thus have agreement with (5.29|) provided we can show 


8 

w L 


dt = 0. 


a=0 


(5.31) 
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Now, using the integral representation from S) shows 

poo -| pOO pi _ -j pOO pi -j 

/ dt=- dt ds Jo{V^)Jq(Vy^) = - dtJo(Vi) ds-M^Ys), 

Jo “=o 4 Jo 4 Jq Jo s 

where the second equality follows by changing variables 1t/s. But 

pOO pOO pOO ^ 

/ Jo(Vi)dt = 2 / sJo(s)ds = 2 / — (s Ji(s)) ds = 0 

Jo Jo Jo ds 


so indeed ( 5.31 ) holds true. 

Let us now return to the consideration of /scaled general. Previous studies have given formulas of 
a different form to (5.27) and ( 5.28| ). These read | p^ 


dt 


fZ^iX,Y) = K'°'^^\X,Y) 

/“ft (X, r) = x-f* (X, r) + i Ai(y) Ai{t) 

i") = ifF) + Ja-l{u) du 


(5.32) 

(5.33) 

(5.34) 


where here /hard(^!^) i® scaled hard edge with singularity as a; ^ 0+. Agreement between 

(^.27) and (^.32|), (^.33[) is immediate upon substituting the integral formulas from (4.11), (4.12) in (5.27) 


and integrating by parts. It remains to show that the RHS of ( 5.34 ) agrees with the RHS of ( 5.28 ) with 
the replacement a i—> a + 1. The verification is done by making use of Lemma in ( ^.28D to obtain 


fLiUx,Y) 


= i^hard/j^ yN-^ 

a^a+1 2 4VF 


1 8 Y 

J„(v^)Ja+l(VF) + -—-X>'-'i(t,F)d< 


1 d 

--—VYJa+l{VY) t-^JaYt)dt 


(5.35) 


(c.f. ( ^.30| )). Using the integral formula from (without the restriction a = 0), a straightforward 

calculation shows 

^ Y) dt = Y) + t-^Ja{Vi) dt. 


Substituting this in ( |5.35D gives agreement with ( ^.34 ) provided 

^ =J,(VX)J,+ l(VF) + l^^i^ J t-lja(Vt)d<- i^VFja+l(VF)^ t-^Ja{Vi)dt 

Ja+l{VY) 


aVy" 


pOO 

■ / Ja-l{u)du. 

Jvx 


4\/F dvX 

Since both sides vanish as X ^ oo, it suffices to check that the derivative with respect to X of both sides 
agrees. This is easily verified using suitable Bessel function identities. 


5.2 The symplectic symmetry limit 


We know from ( |3.16D and ( |3.17| ) that in the limit A ^ — oo the parameter dependent Laguerre and 
Jacobi ensembles tend to the LSE with weight and the JSE with weight (1 — respectively. 


The results of the present study give that the corresponding fc-point distribution is given by (1.31) with 
matrix elements specified by (3.99)-(^.101). As with the orthogonal symmetry limit, previous studies 
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ft hj, y have obtained the fc-point distribution for the LSE with general a > — 1, and the JSE with 


general a,b > —1. In particular, in terms of (5.1) and (p.2D, the results of ly give that the fc-point 
distribution for the matrix ensembles SE„(a::“e“'^) (LSE) and SE„((1 +a;)^(l — a;)“+^) (JSE), as specified 
by (HI), are given by (HI) with the 2x2 matrix / in (|l.3l|) having as its top left entry 

' ^2n,a—1 Jx 

(5.36) 

in the Laguerre case, and 

2/ {x,y) = (l-a;)(Y—) K^n,a,b-i{x,y) - 


-(1 - y ) 


(a-l)/2 


^_i,a,b-i(2’T' +\{a + b-2)) 

J X 

in the Jacobi case (we have taken the transpose of the qdet formula in [y, and thus have interchanged x 
and y in relative to the expression in ||]). The other entries in the matrix / are related to by 

f^{x,y) = f^{y,x), f^{x,y) = -J f^{y,t)dt, f^{x,y) = -^f^{y,x). (5.38) 

We see from the first equality in ( 3.65| ), ( $.100 ) and ( 3.101 ) that the equations ( $.38| ) are identical to 
those obtained in Section 3.4 for the symplectic limit in the Laguerre case. Furthermore, in the Jacobi 
case, defining 

f^{x,y) = ^-^f^{x,y) = ^{l-y)KiJx,y) + ^(^il-x)-^-l^ j Ki^{t,x)dt (5.39) 

the first equality in ( 3.65| ) and (3.99)-( $.101 ) together with simple scaling invariances of the quaternion 
determinant imply that the A:-point correlation can be written as qdet[/] with the elements of the 2x2 
matrix / given as in ( 5.38 ) but with each J'*® replaced by . Thus it remains to show agreement 
between the first equality in ( 3.99| ), substituted in the first equation of ( $.38 ), with ( $.36| ) in the case 
a = 0, and ( $.39 ), substituted in the first equation of ( $.38| ), with ( 5.381) in the case 6 = 0. 

Consider first the Laguerre case. We know from ® 0 that ( |5.36| ) can be rewritten to read 


where 


G[{y) = 


26^ 

^"- 271 , 0-1 


2/ {x,y)=K^^ {x,y)+G^{y)G^{x) 


ya/2-lg-,/2p.-l(^) 


(5.40) 


G{{x) = hx'^Gf, p (2n + a-1) / dt\ 

^ J X ' 

= t°-G-^e-^G{2nLlJ^x) - (2n + a)L5„_i(x)) . 

For the RHS of ( $.99 ) to equal the RHS of ( 5.40 ) substituted in the first equality of ( $.3^ ) in the case 
a = 0, by setting j/ = 0 in both expressions we see that a necessary condition is that 

1 Id 

-K^n{x,0) + 2^/0 

which is obtained by setting j/ = 0 in both expressions. This is just the previously established identity 
( $.88 ). Knowing that both sides agree at y = 0, to show they are equal for all y it suffices to show 
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that their partial derivatives with respect to y agree. In fact equating the partial derivatives gives the 
previously established identity ( 5.13| ). 

To show agreement in the Jacobi case, we require the analogue of (5.40|). 


Lemma 21. The formula (5.31) has the alternative form 

2/“(a;, ?/) = (1 - y) + Gf (y)G^ (x) 

where 

1 


(5.41) 


Gi{y) = 


An 


in + a + b-l q„_i . 




g2{v) = 

+ (& + 2n - 1) dt (1 - 

j (1 - t)G-^^G(l + f'jb/ 2 p{a,b+l) 


a + b + 2n 


Proof. This is derived in a similar fashion to the result of Lemma 18 


□ 

With (5.41) substituted in the first equality of ( 5.38| ), agreement of the resulting expression with ( ^.39 ) 
in the case b = 0 can be established by first checking that both expressions coincide at y = — 1 (this 
follows from ( ^.69 )), and then showing both expressions have the same partial derivative with respect to 
y (this follows from (5.15)). 

Next we consider the form of the scaled matrix elements of Propositions [l^ and ^ as a ^ — oo, which 
corresponds to the symplectic limit. Integration by parts shows 


Lf 

J scaled 




fsciUX^Y) = 


1 Y 


2 \ X 


K 


scaled 


iX.Y) 


1 d 
Ju dX 


^scaled (5.42) 


-f!c\UX,Y) 


J X 

A 

W 


fscb^ledit^Y) dt 


f!clUX,Y) 


in the case of the soft edge and the bulk, and 


Y 


^fii^iX,Y)r^ fii^iX,Y) = ^ 


1 Y 


^hard(^, Y') p 


1 d 
Ju dX 


VXK'^‘^’^^{X,t)dt 


(5.43) 


8 

^XYfZUX, Y)^- fZdit, Y) dt 

a Jo 

{X,Y) 


8 ^ 


in the case of the hard edge with singularity x°‘^^ as a; ^ 0. 

Previous studies have given formulas of a different form to (5.42|) and (^.43|), these being |17 


&iX,Y) = K'^'^^\X,Y) 
f!,l{X,Y) = iL-f*(X,r)--Ai(X)y^ Ai{s)ds 

f^l^iX, Y) = Y) - ^ J,_i(v/J) J„+i(s) ds, 


40 






























where here /hard relates to the scaled hard edge with singularity as a; ^ 0+. The verification that 
these formulas are identical to (5.42) or ( ^.43 ) is done in the same way as reconciling ( ^.27 )-(5.28) with 
(|^ and dJI). 


5.3 Bulk structure function for the decimated orthogonal ensemble with a 
parameter 


In general the 2-point correlation implied by (1.31) is given in terms of the corresponding 1-point corre¬ 
lation (the density) and the matrix elements (s, s' = 1, 2) of / according to 

P2{,Xi,X2) = Pl{xi)pi{X2) - (^f^^{Xi,X2)f'^{xi,X2) - f^'^{xi,X2)f^{xi,X2)^. (5.44) 

Now the truncated quantity 

P^{xi,X2) := P2{xi,X2) - Pi{xi)pi{x2) (5.45) 

decays for large \xi — X 2 \- Furthermore, in the bulk it is a function only of the difference \xi — X 2 \- Thus 
in the bulk the Fourier transform of p^{x\ — X2,0) is a well defined quantity. We seek its evaluation in 
the case of the parameter dependent orthogonal ensemble with a parameter, when the matrix elements 


are specified by (|4.27| )- ([4.29|) . Explicitly, we will compute the quantity 

/ CXD 

pl{x)e^’^^dx 

-OO 

/ CXD 

{f^\x)r{-x) + f^^{x)r{-x))e 

-OO 


dx 


(5.46) 


where to obtain the second equality we have substituted ( 5.44 ) and made use of the facts, apparant from 
(p^-dO^), that /“ (cc), P^{x) are even functions of x, while are odd functions of x. Following 

[p^ we know that the most efficient way to compute such a Fourier transform is to make use of the 
general formula 

/ OO -| poo nOO 

f{x)f{-x)e^>^^dx=— mf{l-k)dl, f{x)e^‘^dx. (5.47) 

-OO J — oo J — oo 

Thus we must first compute the Fourier transform of the individual matrix elements. 

The calculation of the Fourier transform of the individual matrix elements is simplified by Hrst noting. 


making use of the fact that Y) is a function of X — Y, that the expression ( 4.27 ) for fPpX, Y) 

can be simplihed to read 

PIPX,Y)=K^'^^^{X,Y), 

independent of the parameter a. Recalling the first equality in (|3.65 ), and the integral representation of 
( |4.11 ), we thus have 

/bulk(0 = /bulk(0 = X|i|<7r- 

Furthermore, substituting the integral representation of (|4.11 ) in (4.28)-(4.29) allows us to compute 


/bulk(0 = 


(a/2)2 + 

2il 


Thus 


PU{i)KUii-k)di = 


/ P^il)P^il-k)dl = - 

/ —OO J — 7T—\k\ 


$21 

X.\l\<TTi /bulkw/ (o;/2)^ + 

J 27r- |fc|, |fc| < 27r 

1^1 ^2, 

min(— ) (a/2)2 + (Z + |A:|)2 ; 


l+\k\ 


(a/2)2+P 


(5.48) 
dl. (5.49) 
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The integral (Pi can be evaluated in terms of elementary functions, with there being a different 
functional form for |fc| < 27r, |fc| > 27r. Substituting this and (5.48) in (^.46[) (appropriately rewritten 


using (5.47)), we hnd that for |fc| < 27r 


S{k) = + 

TT 


\k\ 


47r((a/2)2 + k^) . 


27r 2|fc| + 

arctan-arctan- 

a a 


,a 




log ( 


ia/2y 


Hfc|( 

+ i\k\+nr 


27T 


) 


,/2y 


a 


- y log|l - |fc|/7r| 


(5.50) 


while for |fc| > 27r 


S{k) = 2 + 


arctan(2(|fc| — 7r)/a) arctan(2(|fc| +7r)/a) 


47r((a/2)2 + k"^) 


47r((a/2)2 + k^) 


(2(a/2)^|fc| + |fc|3) (a/2)^ + (|fc|+^)^ 

4^((a/2)2 + P) °S(a/2)2 + (|fc|_^)2- 


Of particular interest is the small |A:| expansion of (5.50). We find 

^2 ^2 


Sik) = 


\k\ , 1 (a/2)2-^2 ((^/2)2_^2)2 


TT 27r2 (q!/ 2)2 + TT^ 47r3((a/2)2 + Tr^)^ 


|fc|3 + 0(fc4). 


(5.51) 


(5.52) 


We see from (5.52) that the coefficient of the leading order term, proportional to |fc|, in the small 
|A:| expansion of S{k) is independent of the parameter a, and has the value I/tt. This is to anticipated 
from the interpretation of (O) as a one-component log-potential Coulomb system with coupling /3 = 1. 
The coupling within pairs can be regarded as a short range potential which should not affect properties 
determined by the long-ranged logarithmic potential. One such property is the behaviour 


5(fc) 


\k\ 


Ifcl^o 7r/3 

for a one-component log-potential system with coupling /3 (see e.g. 


1), thus implying the leading 


behaviour seen in (5.52). 


5.4 Distribution of odd labelled coordinates for a special parameter 


In the Introduction, attention was drawn to the properties (0), S) of the parameter dependent 
ensembles ( |1.4| ), ( ]1.16 ) relating to the distribution of the even labelled coordinates. Similarly, we noted 
the properties (|l.l(]|), (|1.19|) of the even labelled coordinates in the superimposed parameter dependent 


ensembles (1^), ( 1.17 ). For the special value of the parameter for which the one body factor for the 
even labelled coordinates reduces to a constant, it turns out that the odd labelled coordinates also have 
a distribution which coincides with that of other matrix ensembles. This follows from the following 
integration formulas. 


Lemma 22. Let xi,X 2 , ■ ■ ■, X 2 n be ordered as in and label this ordering X. We have 


lx 


dX2dX4 ' * * dX2n ^k') — \| ^ 2_7 — 1 — 1 ^2k—l^ 


l<j<k<2n 


and 


lx 


dX2dXi- ■ ■ dX2n {.X 2 j - X 2 k) = —W_X2j-l {X2j-1 - X2k-l)- 


(5.53) 


(5.54) 


l<j<k<n 


j—1 
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Proof. To derive (5.53), we use the Vandermonde determinant formula to write 


(xj — Xk) = det[x 

l<j<k<2n 


k—1 ] 

2n+l-j\j,k=l,...,2n- 


The method of integration over alternative variables gives 


J dx 2 dxi ■ ■ ■ dx 2 n [xj — Xk) = det 


l<j''<fe<2n 


\^2n+l-2i 


{2n)\ 


X 2 j-i det 


j=i 


x’^-^ 

X2n+l-2j 

'^X2n+l-2j 


3 = 1,....n 
k = l,...,2n 


2n+l-2j 

(2^)! 


J = 1.n 

k = l,...,2n 


n*2j-idet 


x'^-^ 

^2n-\-l-2j 
{k — ^)^2n-\-l-2j 


k = l,...,2n 


This final determinant is well known to be equal to the product of differences to the fourth power, and 
thus ( |3.53| ) follows. 

A similar, even simpler, computation gives (5.54). □ 

We remark that dH and (HI are equivalent to the special case of ( |1.12[) and ( [1.14| ) in which (/, g) is 
given by the Jacobi weight in ( 1.13 ) with a = 1. 

For the ensembles (|1.4|), (|1.16), it follows immediately from (^.53|) that 


odd(OE2n(/o,/e)) = SE„(/i) 



^g-a;/2gAa;/2 g-x/2g-^a;/2 ^2^-x'^ 


A=-l 


{x > 0) 


A=l-c 


( — 1 < X < 1) 


(5.55) 


(5.56) 


((1 — x)(“ (1 — (1 + x)^(l — x)“ ^) 

in the Laguerre and Jacobi cases respectively. Similarly, for the ensembles (^]^), (|1.17 ), it follows imme- 

(5.57) 


diately from ( 5.54 ) that 

odd(OE„(/o, /e) U OE„(/o, /e)) = UE„(/l) 
where /□, /e are as in ( ^.56 ), while 


h = 


xe Laguerre case 

(1 + x)(l — x)““^, Jacobi case 


(5.58) 


A consequence of ( ^.55 ) and ( 5.57 ) is that the fc-point odd-odd correlation for the ensemble on the 
LHS must coincide with the fc-point correlation for the ensemble on the RHS. This is simple to explicitly 
verify for the relation ( |5.57| ). Then the fc-point distribution on the RHS is given by 

det[A:,^^i(xj,Xi)]yz=i,...,fe, det[(l i(xj,xi)]j-i=i....,fe (5.59) 

in the Laguerre and Jacobi cases respectively. On the LHS the odd-odd /c-point correlation is given by 


det[K^^(xj,xi) h,i=i.....k, det[(l-Xj)A:/„(xj,Xi)]yi=i,....fc, 


A=-l 


(5.60) 


where are given by (2.15) and ( ^.16 ) respectively. Using appropriate Laguerre and Jacobi 

polynomial formulas to evaluate 

J Li{x)dx, ■^e~'"Li{x), J Pi°'’°\x)dx, ^(1 - x)“p/“’°^(x) 

it is readily seen that (^3.60 ) can be reduced to ( 5.59 ). 
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5.5 Gap probabilities and eigenvalue distributions 


Observable quantities for eigenvalue sequences closely related to correlation functions are gap probabilities 
and distribution functions of individual eigenvalues. The gap probability specifies the probability — to 
be denoted E{p\I\ ME) — that an interval I (the gap) in a given matrix ensemble ME contains precisely 
p eigenvalues. With the eigenvalues ordered as in (O), we may choose to observe eigenvalues of a definite 
parity, even or odd labelled. In this case we denote the gap probability as E^''> {p; I;ME), (•) = (e)ven, 
(o)dd, where p now refers to eigenvalues of parity (•) only. For I = (s,oo) with s inside the support of 
ME, it follows from the ordering dLl) that 

(s, oo); ME) = ^i5(2p — 1; (s, oo); ME) + if(2p; (s, oo); ME)^ 

E(®)(p;(s,oo);ME) = (e;(2p; (s, oo); ME) + E;(2p + 1; (s, oo); ME)) (5.1) 

where £^(—1;/; ME) := 0. Consequently we can express (s,oo);ME) in terms of (s,oo);ME)}, 


p 


p-i 


E(2p;(s,oo);ME) = E;(°)(j; (s, oo); ME) - (s, oo); ME) 

j^O j^O 

p p 

E(2p+l;(s,oo);ME) = ^ (j; (s, oo); ME) - ^ E;(°)(j; (s, oo); ME). (5.2) 

j=o j=o 


We will denote the PDF for the distribution function of the kth eigenvalue Xk (with the ordering (1.5)) 
by p{k — l;s;ME) (here fc — 1 is the number of eigenvalues greater than Xk)- A standard formula (see 
e.g. Q) gives that 


p{k — 1; s; ME) = —E(k — 1; (s, oo); ME) + p{k — 2; s; ME), fc > 1 (5-3) 

as 

wherep(—l; s; ME) := 0, so knowledge of {E{p-, (s, oo); ME)}p=o,....fc-i suffices to compute p(k—l; s; ME). 

Let us consider first the gap probabilities and eigenvalue distributions for ME = (LOE„ U LOE„)^ 
and ME = (JOE„ U JOE„)"^ (i.e. the PDFs ( pT^ ) and ( 1.16| )). The identities ( 1.12| ) and ( 1.18| ) tell us that 


E(^^(p; (s, oo); (L0E„ U L0E„)^) 
E(<^^(p; (s,1);(JOE„U J0E„)^) 
p(2k - 1; s; (LOE„ U LOE„)"^) 
p(2k - 1; s; (JOE„ U JOE„)^) 


E(p; (s, oo); LUE„|a=o), 
A(p; (s, 1 ); JUE„|b=o) 
p{k- l;s;LUE„|a=o), 

p{k 1, s, JUEyjl^— q). 


(5.4) 


which we can check are consistent with ( |5.3| ). Note in particular that each quantity in (5.4) is independent 
of the parameter A, and we remark too that each has a known Painleve transcendent evaluation ^ 
il- To specify E^°'> and thus p(2fc; s; ME), E for these matrix ensembles we make use of the standard 
formula relating the gap probability to the correlation functions. 


E^°\p]I]ME) = + J^dxi ■ ■ ■ J^dxk p^^^ixu ■. ■ (5.5) 


We know from (1.24) that is the k x k determinant det[Koo(xi, In this circumstance 

the expression in brackets is just the expansion of the Fredholm determinant of the integral operator 
Koo with kernel Koo(x,y) supported on J, 


^ ^ f-l)P dP 

(p; /; ME) = — det (1 - ^i^oo) [ 5=1 


p\ 


(5.6) 
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(here 1 denotes the identity operator). When p = 0 this reads 


S(°)(0;/;ME) =£;(0;/;ME) = det(l - Koo) 


(5.7) 


where the first equality follows from the first equation of ( ^.l[ ). Using the fact that with A = 0, (LOE„ U 
LOE„)"^ and (JOE„ U JOE„)^ reduce to LOE„ U LOE„ and JOE„ U JOE„ respectively, it follows that 
if(0;/; ME|^=o) = (£'(0;/; OE))^ and thus we deduce from (5.7) that 

(e(0; (s,oo);LOEU=o))" = det(l - K^,\a=o), (^0; (s, oo); JOE| a„(a^-i )/2 )) " = det(l - KHa^o)- 

(5.8) 

We remark that if(0; (s, oo); LOE|a=o) has recently been evaluated in terms of Painleve transcendents 

iil- 

Let us now consider the gap probabilities and eigenvalue distributions for ME = (LOE 2 „)'^ and 
ME = (J0E2„)"^ (i.e. the PDFs (IM) and ( 1.16 )). The identities ( 1.14 ) and ( 1.19 ) tell us that 


B(^)(p; (s,oo);(LOE2„)^) 
(s,1);(JOE2„)^) 
p(2k-l;s; (L 0 E 2 „)^) 
p(2fc- l;s;(JOE„U JOE„)^) 


B(p; (s,oo);LSE„la=o) 

E{p- (s,l);JSE„U„.+i) 
6=0 

p{k - 1; s; LSE„|o=o), 

p(k - 1; s; JSE„| a„a+i), 
6=0 


(5.9) 


(c.f. (^)). Of these quantities ^(0; (s, cx)); LSE^ja^o) and p(0; 5 ; LSE^ila^o) are known in terms of 
Painleve transcendents |^, To specify we again make use of (^), this time noting from ( 1.27 ) 
that is the k x k quaternion determinant qdet[fooixi, where foo is the 2x2 matrix 

representation of a particular real quanternion, which in turn implies 


^ ^ 1 - 1 )^’ dP 

E^°\p ■ J;ME) = qdet(l -efoo)lc=i 

p! 

(with {Aj} denoting the distinct eigenvalues of /oo, qdet(l — C/oo) = 0^(1 ~ We remark that for 

general A, E^°\0; (s,oo), (LOE)^„) has recently been evaluated in terms of Painleve transcendents 
All the above formulas have well defined scaled limits. In particular 


lim E(°\p; (4n + 2(2ny/^s,oo); (LOE„ U LOE„)^=“/ 2 ( 2 n) 


1/3, 


:= e(°^ -f‘(p; (s, oo); (OE U OE)“) = det(l - 




lim E(°)(p; (1 - 1); (JOE„ U JOE„)^=4n=“) 

n—^oo Zn 

f—1)P f)P 

:= ( 0 , s). (OEU OE)“’“) = det(l - 

(for the justihcation of the limiting processes see 1^). An evaluation of (s,oo); (OE U OE)“) 

in terms of a Riemann-Hilbert problem and Painleve II transcendents has been given in . 
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